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in  many  different  systems.  The  propagation  of  azimuthal  spokes  are  investigated  in  a 
6  kW  HET  known  as  the  H6  using  ultra-fast  imaging  and  azimuthally  spaced  probes.  A 
spoke  surface  is  a  2-D  plot  of  azimuthal  light  intensity  evolution  over  time  calculated  from 
87,500  frames/s  videos.  The  spoke  velocity  has  been  determined  using  three  methods  with 
similar  results:  manual  fitting  of  diagonal  lines  on  the  spoke  surface,  linear  cross-correlation 
between  azimuthal  locations  and  an  approximated  dispersion  relation.  The  spoke  velocity 
for  three  discharge  voltages  (300,  400  and  450  V)  and  three  anode  mass  flow  rates  (14.7, 

19.5  and  25.2  mg/s)  yielded  spoke  velocities  between  1500  and  2200  m/s  across  a  range 
of  normalized  magnetic  field  settings.  The  spoke  velocity  was  inversely  dependent  on 
magnetic  field  strength  for  low  B-field  settings  and  asymptoted  at  B-field  higher  values. 

The  velocities  and  frequencies  are  compared  to  standard  drifts  and  plasma  waves  such 
as  E  x  B  drift,  electrostatic  ion  cyclotron,  magnetosonic  and  various  drift  waves.  The 
empirically  approximated  dispersion  relation  yielded  a  characteristic  velocity  that  matched 
the  ion  acoustic  speed  for  ~  5  eV  electrons  that  exist  in  the  near-anode  and  near-field  plume 
regions  of  the  discharge  channel  based  on  internal  measurements.  Thruster  performance 
has  been  linked  to  operating  mode  where  thrust-to-power  is  maximized  when  azimuthal 
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=  amplitude  for  Lorentzian  fit,  arb.  units  Hz-1 
=  magnetic  field,  T 
=  normalized  magnetic  field 
=  speed  of  light,  ms-1 
=  linear  fit  coefficients 
=  electric  field,  V  m-1 
=  frequency,  Hz 

=  camera  frame  rate,  frames  s-1 
=  peak  frequency  for  spoke  order  m,  Hz 
=  center  frequency  for  Lorentzian  fit,  Hz 
=  inner,  outer  magnet  coil  current,  A 
=  discharge  current  density,  A  cm-1 
=  azimuthal  wave  number,  rad  m-1 
=  discharge  channel  length,  m 
=  probe  spacing,  m 

=  gradient  length  scale  for  parameter  p,  m 
=  spoke  order 
=  minimum  spoke  order 
=  ion,  electron  mass,  kg 
=  neutral  density,  m-3 
=  number  of  bins 
=  number  of  frames 
=  plasma  density,  m-3 
=  ion,  electron  density,  m-3 
=  power  spectral  density,  arb.  units  Hz-1 
=  elementary  charge,  C 
=  mean  discharge  channel  radius,  m 
=  radial  location,  m 
=  electron  temperature,  eV 
=  probe  time  delay,  s 
=  time  delay  from  bin  n  to  m,  s 
=  plasma  potential,  V 
=  Alfven  speed,  ms-1 
=  characteristic  velocity,  ms-1 
=  group  velocity,  ms-1 
=  phase  velocity,  ms-1 
=  ion  acoustic  velocity,  ms-1 
=  spoke  velocity,  ms-1 
=  spoke  velocity  from  bin  n  to  m,  m  s_1 
=  azimuthal  velocity,  ms-1 
=  weighting  for  spoke  order  m 
=  axial  location,  m 

=  dispersion  relation  power  dependence 
=  spoke  velocity  to  Br/B*  power  dependence 
=  full- width  at  half  maximum  for  Lorentzian  fit,  Hz 
=  permittivity  of  free  space,  F  m_1 
=  plasma  resistivity,  Q  m 
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A 6j:k  =  angular  difference  from  bin  n  to  m,  deg 

0sp  =  spoke  angular  velocity,  deg  s_1 

A  =  wavelength,  m 

/io  =  permeability  of  free  space,  H  m-1 

crx  =  standard  deviation  or  standard  error  of  variable  x,  units  of  x 

rs  =  shutter  period,  s 

(f)pr  =  probe  phase  lag,  rad 

De  =  electron  Hall  parameter,  rad  s-1 

c J  =  frequency,  rad  s_1 

Uch  =  characteristic  frequency,  rad  s_1 

ooci,e  =  ion,  electron  cyclotron  frequency,  rad  s-1 

c JJih  =  lower  hybrid  frequency,  rad  s_1 

uj pi^  =  ion,  electron  plasma  frequency,  rad  s-1 

I.  Introduction 

Spokes  were  first  observed  in  HET  type  devices  by  Janes  and  Lowder1  and  they  were  identified  as  a 
possible  mechanism  for  cross-field  transport.  Research  on  a  Cylindrical  Hall  Thruster  (CHT)2,3  showed  that 
up  to  half  of  the  discharge  current  can  pass  through  a  spoke.  It  was  also  reported  that  the  CHT  performance 
(measured  by  discharge  current)  increased  when  the  spoke  was  not  present.4  This  can  lead  one  to  conclude 
that  spokes  are  detrimental  to  HET  performance.  However,  the  annular  device  studied  by  Janes  and  Lowder 
differed  from  modern  HETs  and  the  CHT  has  a  significantly  different  magnetic  field  topology  and  physical 
geometry  (no  inner  wall);  so  direct  comparisons  may  not  be  appropriate  with  those  systems  to  modern 
annular  HETs  such  as  the  SPT-100  and  H6. 

Research  by  Brown5  and  McDonald6  on  low  voltage  operation  showed  that  thruster  performance  was 
increased  when  spokes  were  stronger  in  the  H6.  Recent  results  by  Sekerak7  on  mode  transitions  clearly 
shows  that  spoke  behavior  was  dominant  in  so-called  ”  local  oscillation  mode”  where  the  thruster  exhibited 
lower  mean  discharge  current  and  discharge  current  oscillation  amplitude.  The  H6  thrust-to-power  are 
maximum  when  the  thruster  is  operating  in  local  mode  with  spokes  clearly  propagating  and  no  significant 
breathing  mode.  Ref.  7  raises  the  causality  question  of  whether  spokes  are  responsible  for  the  improved 
thruster  performance  or  are  indicators  when  the  thruster  is  running  optimally.  Regardless,  the  association  of 
spokes  with  improved  thruster  performance  in  the  H6  drives  us  to  understand  the  fundamental  mechanisms 
of  spoke  mechanics  in  order  to  improve  thruster  operation. 

II.  Mode  Transition  Oscillations 

HETs  have  been  under  development  for  over  50  years  with  significant  experimental  and  flight  histories8 
and  mode  transitions  have  been  commonly  observed  throughout  their  development  as  noted  by  some  of  the 
early  pioneering  Russian  research.9  HETs  have  several  parameters  that  define  a  single  operating  point  such 
as  discharge  voltage,  magnetic  field  strength  (or  magnet  coil  current),  anode  mass  flow  rate  and  cathode 
mass  flow  rate.  Laboratory  HET  discharge  power  supplies  operate  in  voltage  regulated  mode  where  the 
discharge  voltage  between  the  anode  and  cathode  is  held  constant  and  the  discharge  current  is  allowed  to 
fluctuate.  A  general,  qualitative  description  of  mode  transition  can  be  deduced  from  previous  research  as 
the  point  while  varying  one  parameter  and  maintaining  all  others  constant,  a  sharp  discontinuity  is  observed 
in  the  discharge  current  mean  and  oscillation  amplitude.  In  one  mode,  the  discharge  current  oscillation 
amplitude  is  small  with  respect  to  the  mean  discharge  current  value,  while  after  the  mode  transition  the 
mean  discharge  current  rises  sharply  as  well  as  the  oscillation  amplitude.  Previous  researchers  have  identified 
mode  transitions  in  HETs  where  a  small  change  in  a  thruster  operating  parameter  such  as  discharge  voltage, 
magnetic  field  or  mass  flow  rates  causes  the  thruster  discharge  current  mean  and  oscillation  amplitude  to 
increase  significantly.  Mode  transitions  in  the  H6  were  induced  by  varying  the  magnetic  field  intensity  while 
holding  all  other  operating  parameters  constant  and  measurements  were  acquired  with  ultra-fast  imaging 
and  ion  saturation  probes. 
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A.  Thruster  and  Facilities 


A  recent  experiment  induced  mode  transitions  in  a  6  kW  class  laboratory  HET  called  the  H6  shown  in  Fig.  1. 
The  experimental  setup  is  described  in  detail  in  Ref.  7  and  only  pertinent  details  will  be  recapitulated  here. 
The  investigation  was  conducted  in  the  Large  Vacuum  Test  Facility  (LVTF)  of  the  Plasmadynamics  and 
Electric  Propulsion  Laboratory  (PEPL)  at  the  University  of  Michigan.  The  test  matrix  included  variations 
in  discharge  voltage  and  xenon  propellant  flow  rates.  Propellant  mass  flow  rates  tested  were  25.2,  19.5 
and  14.7  mg/s  through  the  anode  and  1.8,  1.4  and  1.0  mg/s  (7%  cathode  flow  fraction)  through  the  LaE>6 
cathode.  Discharge  voltages  of  300,  400  and  450  V  were  applied  between  the  anode  and  cathode.  The  xenon 
corrected  chamber  pressures  were  6.1  x  10'  6,  9.0  x  10  6  and  1.2  x  10  5  Torr  for  14.7,  19.5  and  25.2  mg/s 
anode  flow  rate,  respectively. 

The  H6  was  a  joint  development  effort  of  the  University  of  Michigan,  the  Air  Force  Research  Laboratory 
(AFRL)  at  Edwards  AFB,  and  the  NASA  Jet  Propulsion  Laboratory  (JPL);  a  separate  copy  of  the  thruster 
is  maintained  at  each  institution.  It  is  notable  for  its  high  total  efficiency;  e.g.,  64%  at  300- V  (6-kW)  with 
a  specific  impulse  of  1950  s  and  70%  at  800-V  (6-kW)  with  a  specific  impulse  of  3170  s.10  Since  first  firing 
in  2006,  the  H6  has  been  well  characterized  by  a  variety  of  diagnostic  and  modeling  techniques  at  Michigan, 
AFRL  and  JPL.  Six  experimental  doctoral  dissertations  have  focused  on  the  thruster  to  date,  which  have 
spawned  numerous  associated  conference  and  journal  articles. 

The  magnetic  field  shape  shown  in  Fig.  1  was 
kept  constant  during  testing,  though  the  magnitude, 
noted  by  Br/B *,  was  varied  throughout  the  testing 
in  order  to  induce  a  mode  change  within  the  H6. 

The  quantity  Br/B*  is  the  maximum  radial  mag¬ 
netic  field  value  at  a  particular  setting  of  inner  mag¬ 
net  (Jim)  current  and  outer  magnet  Jom)  current 
divided  by  the  nominal  maximum  radial  magnetic 
field.  The  nominal  magnetic  field  ( B '*)  strength 
(i.e.,  for  maximum  thruster  performance)  at  300  V 
and  20  A  discharge  current  is  J/m  =  3.50  A  and 
Iom  =  3.13  A.10 

B.  Fast  Cam  and  ISR  Probes 


Figure  1:  (Left)  H6  with  direction  of  magnetic  field 
and  E  x  B  shown.  (Right)  Profile  picture  of  the  H6 
operating  at  nominal  conditions  with  magnetic  field 
stream  lines  overlaid  and  discharge  channel  outlined. 
Reproduced  from  Figure  6  of  Ref.  7. 


High-speed  imaging  was  acquired  with  a  Photron 
SA5  FASTCAM  with  a  Nikon  ED  AF  Nikkor  80-200 
mm  lens  at  its  maximum  aperture  f/2.8.  The  SA5 
is  capable  of  up  to  1,000,000  fps  with  128x16  pixel 
resolution,  but  was  used  at  87,500  fps  with  256x256 
pixel  resolution  for  this  testing.  The  camera  was  6  m 
downstream  from  the  thruster  outside  LVTF  with  a 
view  of  the  thruster  through  a  viewport.  Details  of 
the  McDonald  technique  are  given  in  Refs.  11  and 

7,  but  a  brief  summary  is  provided  here.  The  video  is  imported  into  MathWorks  MatLab  where  each  frame 
is  a  256  x  256  matrix  of  light  intensity  values  and  the  AC  component  is  subtracted  individually  from  each 
pixel.  The  discharge  channel  is  isolated  and  divided  into  180  two-degree  bins.  The  pixels  in  each  bin  are 
averaged  together  generating  a  180  x  1  vector  of  light  intensity  for  each  frame.  A  2-D  plot  of  all  frames  is 
the  spoke  surface  where  the  ordinate  is  azimuthal  location  around  the  discharge  channel  in  clock  positions 
and  the  abscissa  is  time  with  each  vertical  column  of  values  representing  one  frame  of  video.  Adding  all 
bins  together  yields  the  m  =  0  or  mo  mode  and  was  first  shown  by  Lobbia12  to  linearly  correlate  to  the 
discharge  current.  During  the  mode  transition  investigation  of  Ref.  7  the  same  strong,  linear  correlation  was 
also  observed  and  used  as  the  basis  for  converting  light  intensity  to  discharge  current  density  with  several 
assumptions. 

The  spoke  surface  or  discharge  current  density  surface  yields  valuable  information  on  plasma  oscillations 
within  the  discharge  channel  by  showing  the  time-resolved,  azimuthal  distribution  of  light  intensity.  Vertical 
features  represent  extremes  in  discharge  current  density  that  occur  everywhere  in  the  channel  simultane¬ 
ously.  Diagonal  features  are  perturbations  in  discharge  current  density  that  propagate  azimuthally  around 
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the  discharge  channel.  Lines  from  upper-left  to  lower-right  are  propagating  counter-clockwise  around  the 
discharge  channel  and  lines  from  the  lower-left  to  upper-right  are  propagating  clockwise.  The  Ex  B  direction 
in  the  H6  is  counter-clockwise  since  the  B-field  direction  is  radially  out.  It  was  shown  that  all  azimuthally 
propagating  features  are  in  that  direction  represented  by  lines  from  upper-left  to  lower-right  with  the  slope 
corresponding  to  propagation  velocity  in  deg/s.  The  High-speed  Image  Analysis  (HIA)  and  two-dimensional 
Power  Spectral  Density  (PSD)  are  discussed  further  in  Section  III  C. 

Ion  saturation  reference  (ISR)  probes  were  positioned  1.5  discharge  channel  radii  down  stream  at  the  6 
o’clock  position  on  thruster  center  line  a  shown  in  Figure  9  of  Ref.  7.  The  ISR  probe  gap  was  29.5  ±0.5 
mm  apart,  which  corresponds  to  21.4°  ±  1.7°  of  azimuthal  spacing;  i.e.,  ^  11°  on  either  side  of  6  o’clock. 
The  probes  were  0.13- mm  diameter  pure  tungsten  wire  with  3  mm  exposed.  The  ISR  probes  were  biased  to 
-30  V  with  respect  to  ground,  which  is  more  than  16  Te  below  the  plasma  potential  and  therefore  are  safely 
in  ion  saturation.  The  ISR  current  is  measured  external  to  the  LVTF  across  a  100-fl  shunt  resistor  through 
an  Analog  Devices  AD  215  120-kHz  low-distortion  isolation  amplifier.  The  data  acquisition  system  used  to 
record  discharge  current  and  the  ISR  signals  are  sampled  at  rates  up  to  180  MHz  with  16-bit  AlazarTech 
ATS9462  digitizers. 

C.  Transition  Results 

The  magnetic  field  Hr/F>*  was  varied  by  changing  the  inner  and  outer  magnet  coil  currents  in  a  constant 
ratio  with  all  other  parameters  held  constant  including  flow  rates,  discharge  voltage,  and  chamber  pressure. 
Maintaining  a  constant  1.12  ratio  of  inner  to  outer  coil  current  allowed  the  magnetic  field  magnitude  to 
be  varied  without  changing  the  shape  shown  in  Fig.  1.  Decreasing  Br/B*  below  a  certain  threshold  was 
shown  to  repeatedly  induce  a  mode  transition  where  the  mean  discharge  current  increased  and  the  discharge 
current  amplitude  increased.  An  example  of  this  transition  is  shown  at  the  top  of  Fig.  2.  The  mean  discharge 
current  and  oscillation  amplitude  (Root-Mean-Square  or  RMS)  are  lowest  in  local  mode  ( Br/B *  >  0.61  and 
then  increase  sharply  in  global  mode  ( Br/B *  <  0.61).  The  transition  point  is  Br/B*\trans  =  0.61.  A  defined 
transition  Br/B*  is  misleading  because  there  is  a  transition  region  where  the  plasma  exhibits  both  types  of 
oscillations  as  shown  in  Fig.  2,  however  the  transition  typically  occurred  over  only  ~  10%  change  in  Br/ B*. 
Figure  12,  13  and  14  of  Ref.  7  shows  that  Br / B*\trans  increased  with  increasing  flow  rate  and  discharge 
voltage.  The  spokes  shown  in  Fig.  2  for  Hr/F>*  =  1.00,  1.48  are  localized  oscillations  that  are  typically 
10-20%  of  the  mean  discharge  current  density  value,  while  the  oscillations  in  global  mode  Hr/H*  =  0.52  can 
be  100%  of  the  mean  value. 

The  modes  are  described  in  Ref.  7  as  global  oscillation  mode  and  local  oscillation  mode.  In  global 
mode  the  entire  discharge  channel  is  oscillating  in  unison  and  spokes  are  either  absent  or  negligible  with 
discharge  current  oscillation  amplitude  (RMS)  greater  than  10%  of  the  mean  value.  Downstream  azimuthally 
spaced  probes  show  no  signal  delay  between  each  other  and  are  very  well  correlated  to  the  discharge  current 
signal.  In  local  oscillation  mode  perturbations  in  the  discharge  current  density  are  seen  to  propagate  in  the 
E  x  B  direction  with  clear  spokes  shown  in  a  HIA  PSD.  The  discharge  current  oscillation  amplitude  and 
mean  values  are  significantly  lower  than  global  mode.  Downstream  azimuthally  spaced  probes  show  a  clear 
signal  delay  between  each  other  indicating  the  passage  of  spokes  but  are  not  well  correlated  to  the  discharge 
current  indicating  localized  plasma  oscillations  within  the  discharge  channel.  The  mode  transitions  were 
consistent  across  different  tests  and  showed  no  hysteresis,  but  did  change  at  different  operating  conditions. 
The  transition  between  global  mode  and  local  mode  occurred  at  Br/B*  for  higher  mass  flow  rate  or  higher 
discharge  voltage.  The  investigation  did  not  conclude  a  mechanism  that  caused  mode  transitions,  but  that 
work  is  ongoing.  The  thrust  was  constant  within  experimental  error  through  the  mode  transition  but  the 
thrust-to-power  ratio  decreased  by  25%  for  the  14.7-mg/s  flow  rate;  the  peak  in  thrust-to-power  occurred 
near  the  transition  point.  The  plume  showed  significant  differences  between  modes  with  the  global  mode 
significantly  brighter  in  the  channel  and  the  near-field  plasma  as  well  as  exhibiting  a  plasma  spike  on  thruster 
centerline.  For  the  external  cathode  case  the  plasma  spike  disappeared  in  local  oscillation  mode. 

It  was  concluded  that  the  H6  and  likely  any  similar  thruster  should  be  operated  in  local  oscillation 
mode  to  minimize  discharge  current  and  maximize  performance.  Thruster  performance  maps  should  include 
variation  in  discharge  current,  discharge  voltage,  magnetic  field  and  flow  rate  to  identify  transition  regions 
throughout  the  life  of  a  thruster.  Mode  transitions  provide  valuable  insight  to  thruster  operation  and  suggest 
improved  methods  for  thruster  performance  characterization. 
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Figure  2:  B-field  sweep  for  300  V,  19.5  mg/s  showing  transition  at  Br/B*  =  0.61.  The  discharge  current 
mean  and  oscillation  amplitude  are  shown  with  the  transition  and  for  Br/B *  settings  selected  for  further 
analysis.  The  middle  row  plots  are  HIA  PSDs  and  the  bottom  row  plots  are  discharge  current  density.  The 
scale  range  for  Br/B*  =  0.52  discharge  current  density  is  larger  due  to  the  magnitude  of  oscillations.  A 
500-Hz  moving  average  filter  has  been  applied  to  smooth  all  PSDs.  Reproduced  from  Figure  17  of  Ref.  7. 


III.  Spoke  Velocity  Measurements 

The  spoke  velocity  can  be  calculated  from  either  the  high-speed  imaging  or  the  azimuthally  spaced 
probes  using  several  different  methods.  Spokes  are  observed  to  propagate  at  a  range  of  velocities,  so  there 
is  a  distribution  associated  with  the  speed  akin  to  a  distribution  function.  The  methods  below  will  identify 
one  representative  velocity  for  spoke  propagation. 

A.  Manual  Method 

Spokes  are  unambiguously  observed  in  the  FastCam  videos  and  are  obvious  even  to  the  casual  observer  as 
bright  regions  rotating  azimuthally  around  the  discharge  channel.  Using  the  McDonald  technique  to  create 
a  spoke  surface,7, 11  the  spokes  appear  as  diagonal  stripes  in  the  spoke  surface  as  shown  in  Figure  11  of  Ref. 
7.  This  technique  divides  the  discharge  channel  into  180  two-degree  bins  of  averaged  light  intensity  and  a 
video  consisting  of  Nfr  frames  will  yield  a  180  x  Nfr  spokes  surface. 

The  most  obvious  technique  to  calculate  spoke  velocity  is  to  manually  fit  lines  to  the  diagonal  stripes  on 
the  spoke  surface;  the  slope  of  which  represent  spoke  angular  velocity  0sp  in  deg/s.  In  the  FastCam  videos 
and  subsequent  video  enhancement  (c.f.  Figure  2  and  3  from  Ref.  13),  spokes  are  observed  to  fill  the  entire 
channel  width.  Therefore,  spoke  angular  velocity  is  converted  to  a  linear  velocity  using  the  mean  channel 
radius,  Rchni 


\sv  =  (27rRc/mZ/360)  0sp  (1) 

In  order  to  determine  an  average  spoke  angular  velocity,  45  to  50  lines  are  manually  fitted  to  a  normalized 
spoke  surface  as  shown  in  Fig.  3  for  Br/ B*  =  1.00,  300  V  and  19.5  mg/s.  A  normalized  spoke  surface  shows 
the  spokes  more  clearly  without  altering  their  characteristic  slope.  To  normalize  a  spoke  surface,  each 
frame  (vertical  line)  has  its  mean  value  subtracted  and  is  divided  by  its  RMS  value.  In  order  to  test  the 
uncertainty  due  to  human  error  and  repeatability,  50  lines  were  fit  to  the  same  propagating  spoke  with  a 
standard  deviation  of  39  m/s;  this  will  be  shown  to  be  within  the  standard  deviation  of  a  typical  velocity 
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distribution.  The  velocity  distribution  for  the  spoke  surface  example  of  Fig.  3  is  shown  in  Fig.  4  where  the 
spoke  velocity  is  1526  =b  175  m/s  and  the  uncertainty  is  the  standard  deviation  of  the  distribution. 


Figure  3:  One  millisecond  segment  of  a  nor-  Figure  4:  Velocity  distribution  for  the  manu- 

malized  spoke  surface  showing  14  of  47  man-  ally  fitted  lines  in  Fig.  3. 

ually  fitted  lines  for  Br/B *  =  1.00 


More  sophisticated  techniques  will  be  introduced  in  later  sections,  but  those  results  should  be  within 
the  range  of  this  straightforward,  yet  labor-intensive  approach.  Representative  uncertainties  for  the  manual 
method  are  the  mean  uncertainties  in  Fig.  9  of  185  m/s  and  177  m/s  for  300  V  and  400  V,  respectively. 


B.  Correlation  Method 


Starting  with  the  normalized  spoke  surface  as  discussed  above,  the  time-history  signal  of  light  intensity  for 
each  bin  is  a  1  x  Nfr  vector  representing  light  fluctuations  at  that  azimuthal  location  for  the  duration  of  the 
video  which  is  typically  150  to  250  ms  ( Nfr  ^  13  x  103  to  22  x  103).  A  1  ms  segment  of  four  normalized  light 
intensity  traces  are  shown  in  Fig.  5  for  reference.  A  linear  cross-correlation  analysis  (Eq.  10  from  Ref.  7)  of 
the  signals  between  two  bins,  bj  and  bk,  at  different  azimuthal  locations  with  an  angular  difference  of  A 0j^k 
degrees  will  yield  the  time,  tj,k,  it  took  on  average  for  a  spoke  to  propagate  around  the  channel  from  bj  to 
bk.  Fig.  5  shows  an  example  of  the  time  offset  for  three  azimuthal  locations  (30°,  50°  and  70°)  referenced  to 
12  o’clock  on  the  thruster  face  calculated  from  linear  cross-correlation.  Five  peaks  in  light  intensity  (spokes) 
are  selected  and  shown  how  they  propagate  around  the  thruster  in  Fig.  5.  The  spoke  velocity,  vSpj,k,  from 
bj  to  bk  is 

vspj,k  =  (27rRchni/360)A6j,k/tj,k  (2) 

The  spoke  velocity  for  the  correlation  method  is  the  mean  spoke  velocity  calculated  between  Nkins  compared 


v 


sp 


1 

Nbins 


EE  yspj,k 

j  k 


(3) 


In  principle,  the  spoke  velocity  can  be  calculated  from  the  average  time  delay  using  every  combination 
of  the  180  bins,  which  would  be  over  32,000.  However,  practical  considerations  limit  the  range  of  bins  that 
can  be  compared.  The  camera  frame  rate  is  87,500  frames/s  so  each  frame  represents  11.4/is.  A  spoke 
traveling  at  2000  m/s  will  travel  16°  or  8  bins  in  the  time  span  of  one  frame.  Therefore,  a  practical  lower 
limit  is  A0j  ,k  >  20°  or  10  bins.  A  single  spoke  typically  propagates  one-quarter  of  the  discharge  channel 
circumference  for  most  B-field  settings.  In  strong  spoke  regimes,  a  single  spoke  will  propagate  one-half  to 
even  the  entire  channel  circumference.  A  reliable  upper  limit  for  automated  processing  is  to  only  compare 
bins  where  A0jjk  <  70°  or  35  bins.  In  Fig.  5,  6  cycles  can  be  identified  in  ~0.55  ms  which  corresponds 
to  ~11  kHz  or  a  spoke  period  of  rsp  ~  90  fi s.  Due  to  signal  noise,  the  cross-correlation  peak  occasionally 
matches  to  a  spoke  ahead  or  behind  the  correct  spoke  so  the  calculated  offset  time  is  in  error  by  one  or  two 
rsp.  Although  this  occurs  more  often  when  A 6jik  >  90°,  it  occasionally  occurs  for  A 6jjk  <  70°.  These  points 
are  easy  to  identify  via  manual  inspection,  but  reject  criteria  is  set  for  automated  data  processing  so  any 
spoke  velocity  outside  of  500  to  3500  m/s  is  rejected.  In  order  to  reduce  computational  time,  only  90  bins 
(every  other  bin)  are  used  for  reference  start  points.  All  bins  from  bin  10  to  35  CCW  from  the  reference 
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bin  are  used  for  comparison.  Therefore,  j  is  1,  3,  5,  ...  to  180  and  k  is  10  to  35  in  Eq.  (3)  which  yields  a 
maximum  of  Ngins  =  2430  possible  points.  The  spoke  velocities  for  smaller  A 0j^  will  have  larger  uncertainty 
because  half  of  the  camera  frame  period  (5.7 fis)  represents  a  large  fraction  of  the  spoke  travel  time  14/is 
to  travel  20°).  The  standard  deviation  can  be  reduced  by  choosing  a  larger  value  for  the  lower  limit  of  A  6j^ 
instead  of  20°,  but  the  number  of  points  used  in  the  calculation,  will  also  be  reduced,  so  a  balance 

must  be  reached. 

Using  the  correlation  method  on  the  spoke  sur¬ 
face  in  Fig.  3  yields  a  spoke  velocity  of  \sp  = 

1467  zb  272  m/s  where  the  uncertainty  is  the  stan¬ 
dard  deviation  of  the  velocities  used  in  Eq.  (3).  This 
matches  very  well  with  the  manual  technique  de¬ 
scribed  above.  For  this  datum  point  2266  of  the 
possible  2430  points  are  used  for  20°  <  A 0jjk  <  70°. 

The  manual  method  and  correlation  method  pro¬ 
duce  very  similar  results  as  shown  later  in  Figs.  7 
and  9.  Representative  uncertainties  for  the  correla¬ 
tion  method  are  the  mean  uncertainties  in  Fig.  9  of 
277  m/s  and  263  m/s  for  300  V  and  400  V,  respec¬ 
tively.  The  correlation  method  is  important  since  it 
is  an  automated  and  reliable  procedure  of  providing 
the  same  results  as  the  laborious  manual  method. 

C.  Dispersion  Relation  Method 

The  high-speed  image  analysis  (HIA)  method  devel¬ 
oped  by  McDonald14  and  described  in  detail  in  Ref. 

7  generates  Power  Spectral  Densities  (PSD)  from 
the  2-dimensional  spoke  surface.  Fig.  6.  shows  ex¬ 
amples  for  the  300  V,  19.5  mg/s  test  case  where  peaks  are  clearly  visible  for  each  spoke  order,  mn.  As 
described  by  McDonald  in  his  original  derivation,11  m  is  analogous  to  number  of  wave  lengths  per  channel 
circumference.  Hence  m  =  0  or  mo  is  no  wave  in  the  channel  (the  entire  channel  is  dark  or  bright),  m  =  1 
means  one  wave  in  the  channel  (one  half  bright,  the  other  dark),  m  —  2  is  two  waves  per  channel  (two  bright 
regions,  two  dark  regions),  m  =  3  is  three  waves  per  channel  (three  bright  regions,  three  dark  regions),  etc. 
In  the  literature  m  is  often  called  the  wave  mode,  but  we  call  it  spoke  order  to  avoid  nomenclature  confusion 
with  the  HET  operational  modes  discussed  in  Section  II.  The  azimuthal  wave  number,  he,  is  calculated  from 
the  spoke  order  kg  =  m/r.  Fig.  6  shows  that  each  spoke  order  has  a  unique  peak  frequency  that  is  typically 
3-5  kHz  higher  than  the  previous  m.  Therefor,  the  HIA  PSDs  can  be  used  to  generate  dispersion  plots  of 
peak  frequency  u  versus  wave  number  kg. 

The  HIA  PSDs  are  a  powerful  tool  for  understanding  the  plasma  oscillations  associated  with  HET  oper¬ 
ation.  Ref.  15  showed  that  ion  saturation  reference  probes  identified  the  same  peak  frequencies  as  the  HIA 
PSDs  indicating  that  spoke  related  oscillations  extended  out  into  the  plume.  The  same  spoke  surface  that 
generated  the  normalized  spoke  surface  in  Fig.  3  was  used  to  generate  the  HIA  PSD  for  Br/B*  =  1.00  in 
Fig.  6.  Note  the  most  dominant  peak  for  Br/B*  =  1.00  is  m  =  4  at  10.4  kHz,  which  is  close  the  the  crudely 
estimated  frequency  from  Fig.  5  of  ~11  kHz. 

In  order  to  automatically  identify  the  peak  frequencies  for  each  spoke  order,  the  PSDs  are  first  smoothed 
with  a  250  Hz  moving  average  filter  and  the  maximum  value  identified.  Due  to  noise  from  the  DFT  the 
maximum  value  is  not  always  the  frequency  at  the  center  of  the  peaks  seen  in  Fig.  6  so  a  Lorentzian  a  of 
the  form 


,^VvVVvVUv= 


-0  (12  oclock) 

-  30°  CCW,  t  =22.9  |iis 


-50  CCW,/  =  45.7 |is 


Time,  ms 

Figure  5:  Light  intensity  traces  for  4  azimuthal  lo¬ 
cations  from  the  normalized  spoke  surface  in  Fig.  3. 
Selected  locations  are  thruster  12  o’clock  as  the  refer¬ 
ence  and  30°,  50°  and  70°  CCW  from  12  o’clock.  The 
offset  times  are  calculated  via  linear  cross-correlation 
from  12  o’clock  to  the  other  locations.  Five  peaks  have 
been  selected  to  demonstrate  how  spokes  propagate 
CCW  around  the  thruster  using  the  calculated  offset 
times. 


PSD  {f) 


A)  !r 

77  (/-/o)2  +  (ir)2 


(4) 


is  fit  to  a  segment  of  the  PSD  around  maximum  value.  The  fit  variables  in  Eq.  4  are  the  full-width  at  half 
maximum  T,  amplitude  Aq,  and  center  frequency  /q.  This  identifies  the  frequency  /o  or  cjo  at  the  center 
of  the  primary  peak  for  each  m.  Example  dispersion  plots  using  this  technique  for  the  three  HIA  PSDs  are 
also  shown  in  Fig.  6. 


aRecommendation  of  M.  McDonald  for  determining  peaks  in  2-D  FFTs. 
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Frequency,  kHz 


Figure  6:  HIA  PSDs  for  Br/ B*  =  1.25,  1.00  and  0.73  for  300  V  and  19.5  mg/s  with  m  =  0  —  10  shown.  A 
500  Hz  moving  average  window  has  been  applied  to  each  PSD  trace  to  reduce  noise.  The  peak  frequencies 
are  identified  and  plotted  versus  wave  number  for  corresponding  dispersion  relations. 


For  high  magnetic  field  strength,  Br/ B*  =  1.25,  the  higher  spoke  orders  are  most  prominent  with  m  =  10 
showing  a  peak  near  the  same  height  as  m  =  5.  At  the  nominal  setting,  Br/B*  =  1.00,  the  spoke  orders 
m  =  3  —  5  are  an  order  of  magnitude  higher  than  m  >  6,  although  peaks  are  visible  up  to  m  =  10.  At 
the  lowest  magnetic  field  setting,  Br/B*  =  0.73,  spoke  orders  m  =  3  —  5  are  still  dominant  but  lower  in 
magnitude  than  Br/ B*  =  1.00.  Although  very  weak,  peaks  are  still  visible  for  m  =  6  —  9.  As  magnetic  field 
is  increased,  the  frequency  of  each  spoke  order  decreases  as  shown  in  the  PSDs  and  dispersion  plot  in  Fig.  6. 

Dispersion  relations  can  be  of  any  functional  form  and  can  be  quite  complicated  as  will  be  discussed  in 
the  next  section.  However,  we  begin  the  analysis  of  HIA  dispersion  plots  like  the  ones  shown  in  Fig.  6,  with 
a  simple  functional  from 

cja=ci/^  +  c2  (5) 


where  a  is  a  power  dependance  and  c\  and  c2  are  linear  fit  coefficients  from  a  least-squares  curve  fit  to  the 
dispersion  plots.  Here  only  powers  of  a  =  1  and  2  will  be  considered,  but  higher  powers  are  not  precluded. 
In  the  spirit  of  the  many  elementary  plasma  waves  found  in  homogenous  plasmas,  Eq.  (5)  can  be  written  as 
a  dispersion  relation  of  the  form 


L jJ 


ra  ba 
*chK0 


■  V, 


ch 


(6) 


where  vch  =  is  a  characteristic  velocity  such  as  ion  acoustic  speed  or  Alfven  speed  and  uch  =  |c2|1/Q!  is 
a  characteristic  frequency  such  as  the  ion  cyclotron  frequency  or  plasma  frequency.  The  unexpected  minus 
sign  in  Eq.  (6)  results  from  the  fact  that  c2  is  always  less  than  zero  for  all  fits.  This  can  be  seen  from  the 
example  dispersion  plots  in  Fig.  6  where  the  ordinate  intercept  is  cj  <  0  when  extrapolating  backwards  for 
m  <  2  using  the  points  from  3  <  m  <  12.  The  physical  implications  of  c2  <  0  is  a  limit  of  Vchke  >  u ch  for 
( jj  to  be  real,  otherwise  it  will  have  a  growing  imaginary  component  and  thus  be  unstable.  This  implies  the 
only  spoke  orders  that  can  exist  are 

bbl  ^ >  Rchnl  —  Tlflmin  (7) 

^ch 

In  practice,  rnm^n  is  typically  3  or  4.  The  phase  velocity,  v^^,  and  group  velocity,  v^r,  from  the  dispersion 
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relation  in  Eq.  (6)  are 


ot-\  1  /ol 


(8) 


CJ 


Vph  — 


ke 


duj 

V9r  =  dke=  Vph 


Vch\ 

VphJ 


a 


(9) 


Eq.  (8)  shows  that  the  phase  velocity  will  always  be  less  than  the  characteristic  velocity  and  Eq.  ( 9)  shows 
the  group  velocity  will  always  be  greater  than  the  phase  velocity.  In  the  limit  of  (uch/{kdvch))a  1  that 
follows  from  Eq.  (7),  a  binomial  expansion  of  Eq.  (8)  yields  a  simplified  phase  velocity 


Vph  —  Vch 


a 


01 ch 

kgvch 


a 


/  TXlrnin  \  a 

v  m  ) 


(10) 


With  the  FastCam  frame  rate  at  87,500  frames/s  the  Nyquist  limit  is  43.75  kHz  (2.75  x  105  rad/s)  which 
is  the  asymptotic  peak  value  for  m  >  12  observed  in  the  dispersion  plots  of  Fig.  6.  In  fitting  the  simple 
dispersion  relation  of  Eq.  6  to  the  data  in  Fig.  6,  a  parametric  study  was  done  to  determine  the  limits  on 
m.  Three  different  ranges  were  selected  for  spoke  orders:  3  <  m  <  8,  9, 10.  In  general  the  results  were  not 
sensitive  to  the  upper  limit  of  m  used,  but  the  m  =  8  case  had  more  variation  in  characteristic  velocity.  For 
all  future  comparison  plots  the  range  of  m  used  for  curve  fitting  will  be  m  =  3  —  9  and  a  =  1,  2  in  Eq.  (5). 

The  manual  and  correlation  methods  described  above  both  identify  a  single,  dominant  spoke  velocity  for 
a  given  magnetic  field  setting  and  operating  condition.  However,  the  phase  velocity  from  Eq.  (8),  which  is 
assumed  to  be  the  spoke  velocity,  is  a  function  of  wave  number.  Fig.  7  shows  vph  as  a  function  of  spoke 
order  m  as  a  proxy  for  kg  for  300  V  and  400  V  and  a  =  1,  2  (a  =  3  is  very  similar  to  a  =  2  and  is  not 
shown).  A  single,  representative  spoke  velocity  can  be  calculated  from  a  weighted  average  spoke  velocity 
using  the  PSD  value  at  the  peak  frequency  fm  as  the  weighting  factor  wm  for  each  m.  The  spoke  velocity 
and  weighting  factors  are 


9 

vsp  =  E  W™vph 

m= 5 


PSDjfrn ) 

E  PSD(fm) 


(11) 

(12) 


The  HIA  PSDs  shown  in  Fig.  6  show  that  certain  spoke  orders  are  dominant  at  different  magnetic  field 
settings,  with  spoke  orders  m  =  4  and  5  dominant  for  Br/B*  <  1.0.  Fig.  7  shows  for  300  V  the  phase 
velocities  for  m  =  3  are  far  too  low  and  for  400  V  the  the  phase  velocities  for  m  =  3  and  4  are  too  low.  The 
higher  spoke  orders  are  either  dominant  or  the  same  magnitude  as  m  =  4,  5  for  the  higher  magnetic  field 
settings.  The  weighting  methode  of  Eq.  12  accounts  for  the  higher  spoke  order  dominance  at  higher  Br/B* 
values  and  causes  the  upward  shift  above  ^  1,  which  tracks  very  well  with  the  spoke  velocities  calculated  via 
the  manual  and  correlation  method  and  builds  confidence  in  the  dispersion  method.  The  minimum  spoke 
order  m  =  5  was  chosen  in  Eqs.  (11)  and  (12)  such  that  Eq.  (7)  is  satisfied  for  all  conditions.  Fig.  7 
shows  the  velocity  from  the  correlation  method  typically  follows  the  phase  velocity  for  m  =  5  or  6  closely 
for  Br/B*\trans  <  Br/B*  <  1.0  and  for  Br/B*  >  1.0  the  correlation  velocity  follows  m  >  6.  Therefore 
using  the  phase  velocities  for  5  <  m  <  9  to  calculate  a  representative  spoke  velocity  is  reasonable.  Using 
4  <  m  <  12  yields  the  same  shape,  but  shifted  lower  by  100-200  m/s.  Fig.  7  shows  for  300  V  both  a  = 
1  and  2  yield  spoke  velocities  from  the  dispersion  method  that  match  the  manual  and  correlation  methods, 
but  a  =  2  is  better  correlated  for  400  V. 

The  standard  error  aCl  and  aC2  for  fit  coefficients  c\  and  c 2  in  a  linear,  least-squared  fit  are  easily 
calculated  according  to  Ref.  16.  Using  the  error  propagation  equation,17  the  uncertainty  in  phase  velocity 
from  Eq.  (8)  for  wave  number  ke  is 


°’vph 


Vph 

a 


( (El) 2  f IeU  +  ( (El)2  ( 

V  C1  /  \Vph)  V  c2  /  V  yph 


(13) 


The  total  uncertainty  in  spoke  velocity  is  calculated  using  the  weighted  average  technique  of  Eqs.  (11)  and 
(12).  Representative  uncertainties  for  the  dispersion  relation  method  are  the  mean  uncertainties  in  Fig.  9 
of  115  and  47  m/s  for  a  =  1  and  2  at  300  V,  and  157  and  41  m/s  for  a  =  1  and  2  at  400  V. 
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Figure  7:  Comparison  of  phase  velocities  and  spoke  velocities  for  (a,c)  300  V  and  (b,d)  400  V,  19.5  mg/s  for 
a  =  1  (a,b)  and  2  (c,d).  Colored  lines  represent  the  phase  velocity  for  each  spoke  order  m  calculated  with 
Eq.  8.  Red  lines  with  squares  are  spoke  velocities  calculated  with  the  dispersion  method  and  Eq.  11.  Solid 
black  lines  with  circles  are  spoke  velocities  calculated  using  the  correlation  method.  Dashed  black  lines  with 
triangles  are  spoke  velocities  calculated  using  the  manual  method. 


D.  Probe  Dispersion  Relation 

Lobbia  initially  showed  that  light  intensity  oscillations  in  the  discharge  channel  were  linearly  related  to 
electron  density  oscillations  from  3  to  11  discharge  channel  radii  downstream  in  a  BHT-600.12  McDonald 
conducted  research  on  a  segmented  anode  H6  where  frequency  oscillations  in  discharge  current  to  the  seg¬ 
mented  anode  correlated  to  oscillations  observed  with  high-speed  imaging. 11  Recent  research  on  the  H6 
demonstrated  that  down  stream  probes  observe  the  same  frequencies  in  plasma  oscillations  as  the  Fast  Cam 
observes  with  light  intensity  oscillations  in  the  channel.7, 15  Therefore,  plasma  oscillations  in  the  channel 
are  related  to  oscillations  observed  downstream.  In  order  to  justify  using  a  dispersion  relation  for  the  HI  A 
results,  a  dispersion  plot  for  downstream  probes  should  be  related  to  HI  A  dispersion  plots. 

Lobbia  developed  a  frequency  domain  transfer  function  technique  in  order  to  calculate  aggregate  maps  of 
electron  density  oscillations.18  The  same  technique  is  used  to  calculate  a  frequency  domain  transfer  function 
between  the  ISR  probes  that  yields  gain  and  phase  lag  as  a  function  of  frequency.  Note  that  we  have 
employed  set  averaging  as  discussed  by  Lobbia, 19  which  increases  the  signal-to- noise  ratio  and  helps  smooth 
out  turbulent  uncertainties.  Unfortunately,  this  reduces  the  number  of  points  for  the  Fourier  transform  by  the 
number  of  sets  averaged  which  ultimately  decreases  the  frequency  resolution.  For  this  analysis,  the  number 
of  sets  averaged  is  75.  The  local  azimuthal  wave  number  at  the  probes  is  calculated  from  the  phase  lag,  </>pr, 
using  the  probe  gap  Lpr,  by  kg  =  <fipr/Lpr.20  Fig.  8  plots  frequency  versus  azimuthal  wave  number  for  three 
conditions  at  300  V  (same  conditions  shown  in  Fig.  6)  and  one  condition  at  400  V  for  19.5  mg/s  anode  flow 
rate.  Remarkably,  the  dispersion  plots  from  probes  1.5  Rchni  downstream  show  striking  similarities  to  the 
dispersion  plots  calculated  from  HI  A  analysis  described  in  Section  IIIC;  the  correlation  appears  strongest  for 
Br/B*  =  1.25  and  1.00.  Note  there  appears  to  be  an  offset  where  the  probe  dispersion  plots  are  a  slightly 
higher  frequency  by  ~  2  x  104  rad/s  or  ~  3  kHz.  This  indicates  that  the  same  oscillations  observed  in  the 
discharge  channel  with  high-speed  imaging  are  also  detected  in  the  plume  as  shown  by  PSD  analysis  in  Ref. 
15.  This  qualitatively  builds  confidence  in  the  use  of  a  dispersion  relation  to  represent  plasma  oscillations 
from  HI  A. 
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Figure  8:  Probe  dispersion  plots  for  Br/B*  =  1.25,  1.00  and  0.73  at  300  V,  19.5  mg/s  and  Br/B*  —  1.00 
for  400  V,  19.5  mg/s.  The  points  are  the  dispersion  plots  from  the  phase  lag  between  probes  positioned  1.5 
Bchni  downstream  as  a  function  of  frequency  and  the  red  lines  are  the  dispersion  plots  from  HIA  described 
in  Section  IIIC  and  shown  in  Fig.  6. 


E.  Probe  Delay  Method 

ISR  probes  are  used  to  measure  plasma  oscillations  in  the  plume  that  correlate  to  light  intensity  oscillations 
in  the  discharge  channel.  As  discussed  at  length  in  Ref.  7,  both  ISR  probes  observed  the  same  plasma 
oscillations,  but  in  local  mode  the  signal  was  delayed  whereas  in  global  mode  the  oscillations  occurred  nearly 
simultaneously  at  each  probe.  The  time  delay,  td,  was  determined  from  a  linear  cross-correlation  technique 
described  in  Ref  7.  As  shown  in  Figure  20c  of  that  reference  the  time  delay  in  local  mode  was  between  10 
and  15  /is. 

The  spoke  velocity  can  be  calculated  from  the  linear,  azimuthal  distance  between  each  probe,  Lpr  = 
(27rRc/in;/360)A^i?2,  divided  by  the  time  delay 


—  Lpr/^d  (14) 

The  uncertainty  in  Lpr  is  calculated  from  the  probe  spacing  uncertainty  of  1.7°  and  the  uncertainty  of  td  is 
assumed  to  be  10%  of  the  value.12  These  uncertainties  are  used  to  calculate  the  maximum  and  minimum 
values  for  spoke  velocity  at  a  given  setting  vSp\min  =  (Lpr  A  &Lpr)  /  (td  T  atd)-  Representative  uncertainties 
for  the  probe  delay  method  are  the  mean  uncertainties  in  Fig.  9  of  394  m/s  and  419  m/s  for  300  V  and  400 
V,  respectively. 

F.  Spoke  Velocity  Comparison 

Fig.  9  shows  the  spoke  velocity  calculated  via  all  four  methods  discussed  above  with  error  bars  for  300  V  and 
400  V  at  19.5  mg/s  flow  rate.  The  manual,  correlation  and  dispersion  methods  are  all  very  well  correlated. 
The  spoke  velocity  from  probe  delay  is  consistently  higher  by  ^  30%  for  both  conditions  with  the  400  V 
condition  showing  an  unusual  rise  for  Br/B*  >  0.9.  The  reason  for  this  divergence  is  unknown.  The  spoke 
velocity  is  initially  inversely  dependent  on  Br/B*  until  ^  1  until  it  levels  out  for  higher  magnetic  field 
strength.  The  inverse  dependence  of  vsp  on  Br/B*  is  stronger  for  the  300  V  condition  than  400  V. 
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Fig.  10  shows  the  characteristic  velocities  and  mm^n  for  the  dispersion  method  for  300  V,  400  V  and  a  = 
1  and  2.  The  characteristic  velocities  are  higher  for  the  a  =  1  and  for  a  =  2  they  show  the  same  inverse 
dependence  on  Br/B *  until  ^  1,  after  which  the  become  level  at  the  same  value.  The  minimum  spoke  order 
appears  to  be  linearly  dependent  on  Br/B*  with  a  =  2  higher. 

Fig.  11  shows  a  comparison  of  spoke  velocities  calculated  from  the  correlation  method  for  all  five  condi¬ 
tions  tested.  The  300  V,  19.5  mg/s  condition  is  the  average  of  four  sweeps  and  the  400  V,  19.5  mg/s  condition 
is  the  average  of  two  sweeps.  All  conditions  show  the  same  trend  of  spoke  velocity  inversely  dependent  on 
Br/B*  until  ~  1.  Power  dependencies  of  vsp  oc  ( Br/B*)~ 0,5  and  vsp  oc  ( Br/B*)~° ,25  are  shown  in  Fig.  11 
for  reference  purposes  only  and  were  not  generated  from  curve  fits.  The  300  V,  19.5  mg/s  and  14.7  mg/s 
conditions  show  the  strongest  inverse  dependence  closer  to  the  -0.5  reference,  all  others  are  closer  to  the 
-0.25  reference.  For  Br/B *  >  1,  the  300  V,  14.7  mg/s  condition  still  decreases,  but  not  as  steeply  and  300  V, 
19.5  mg/s  actually  increases  velocity  before  stabilizing.  All  other  conditions  are  essentially  constant  for  the 
higher  magnetic  field  settings.  With  the  exception  of  300  V,  14.7  mg/s,  all  conditions  asymptote  between 
1600  and  1700  m/s  for  the  maximum  magnetic  field  settings. 


Figure  9:  Comparison  of  spoke  velocity  calcu¬ 
lation  methods:  manual,  correlation,  disper¬ 
sion  relation  with  a  =  1,2  and  probe  delay 
method  for  (a)  300  V  and  (b)  400  V.  Not  all 
error  bars  are  shown  for  clarity.  For  the  dis¬ 
persion  relations,  m  >  5  has  been  used  in  Eqs. 
(11)  and  (12). 


Figure  10:  Comparison  of  the  (a)  characteris¬ 
tic  velocity  v^  and  (b)  minimum  spoke  order 
from  Eq.  (7)  for  300  V  and  400  V,  19.5  mg/s. 
Power  dependence  a  =  1  and  2  are  considered. 


Combining  the  above  discussion,  we  can  state  the  following  observations  regarding  spoke  velocities;  any 
theory  on  spoke  mechanisms  and  propagation  should  account  for  these  results. 

1.  Predict  a  spoke  velocity  that  is  1500-2200  m/s  in  the  H6;  spoke  velocity  dependence  on  thruster  size  is 
unknown  although  likely  weak  as  noted  by  McDonald.6  Spoke  velocities  are  not  dependent  on  discharge 
voltage  or  mass  flow  rate  to  within  experimental  error.  Propagation  is  in  the  E  x  B  direction. 

2.  The  dispersion  relation  can  be  approximated  by  a  power  law  dependence  uoa  ~  where 

a  >  1.  The  spoke  velocity  vsp  is  less  than  the  characteristic  velocity  vc^  and  is  dependent  on  the 
dominant  spoke  orders,  typically  m  >  4.  In  general,  the  dominant  spoke  order  increase  with  increasing 
magnetic  field  strength. 
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3.  Spoke  velocity  should  be  weakly,  inversely  dependent  on  magnetic  field  for  Brj B*\trans  <  Br/B*  <  1. 
An  example  dependence  of  vsp  oc  B~ &  where  0.25  <  (3  <  0.5  is  shown  in  Fig.  11,  but  other  functional 
forms  are  possible.  For  Br/B*  >  1  the  spoke  velocity  should  nearly  asymptote  to  a  constant  value. 


Figure  11:  Spoke  velocity  calculated  with  the  correlation  method  for  all  conditions  tested.  Parenthetical 
numbers  are  the  number  of  B-field  sweeps  averaged  together.  Reference  lines  for  possible  functional  forms 
of  vsp  dependence  on  Br/ B*  are  shown  for  discussion  purposes  only. 


G.  Spoke  Velocity  Observation  Limits 


Each  spoke  order  represents  the  number  of  light  and  dark  regions  in  the  thruster  so  the  wavelength  is 
A  =  2'KRc}lnilvn-  If  the  spoke  travels  1/2  wavelength  during  the  period  of  time  the  shutter  is  open,  then  the 
bright  region  will  travel  over  the  dark  region  rendering  the  spoke  unobservable  by  the  camera.  Assuming 
the  open  period  of  the  shutter  is  rs  =  1  / fc  where  fc  is  the  camera  frame  rate,  then  the  observable  spoke 
velocity  is 


vsp  < 


ItRchnlfc 

m 


(15) 


For  m  =  3,  5,  7  and  9,  Eq.  15  yields  maximum  observable  spoke  velocities  of  7300,  4400,  3100  and  2400  m/s, 
respectively.  This  is  within  all  of  the  spoke  velocities  for  each  spoke  order  as  shown  in  Fig.  7;  however,  the 
highest  spoke  orders  are  close  to  the  limit. 


IV.  Plasma  Oscillations 

Now  that  we  have  identified  some  important  characteristics  of  spoke  propagation,  we  seek  to  understand 
the  fundamental,  underlying  mechanisms.  We  start  with  the  empirically  identified,  approximate  dispersion 
relation.  The  characteristic  velocity,  vc^  is  shown  in  Fig.  10  to  be  between  2200  and  2800  m/s  for  a  =  1 
and  between  1800  and  2200  for  a  =  2.  In  addition,  Fig.  11  shows  the  spoke  velocity  or  phase  velocity  is 
between  1500  and  2200  m/s  for  all  tested  operating  conditions.  The  characteristic  frequency,  c is  of  the 
same  order  and  usually  a  few  times  larger  than  the  ion  cyclotron  frequency,  wCi.  We  can  now  attempt  to 
identify  any  frequencies  in  the  10’s  kHz  range  and  characteristic  velocities  or  phase  velocities  that  are  in  the 
range  of  ^  2000  m/s.  The  following  analysis  and  discussion  assumes  +1  ions,  but  multiply  charged  ions  are 
known  to  exist  in  HETs.  E  x  B  probe  measurements  in  the  plume  of  the  H6  showed  Xe2+  at  20%  and  Xe3+ 
at  10%  of  Xe1+,  respectively.10 
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A.  Spoke  Mechanisms  and  Locations 

A  coherent  theory  with  supporting  experimental  evidence  for  spoke  propagation  has  yet  to  emerge,21  so 
the  mechanism  and  location  in  the  plasma  for  their  formation  and  propagation  are  unknown.  However, 
there  have  been  no  shortage  of  research  published  on  azimuthal  oscillations  since  the  early  work  of  Morozov, 
Tilinin  and  Espipchuk  in  the  former  Soviet  Union  and  Janes  and  Lowder  in  the  United  States.  Escobar 
recently  provided  an  excellent  overview  of  experimental,  theoretical  and  numerical  research  on  azimuthal 
oscillations,22  and  one  should  reference  that  work  for  a  comprehensive  discussion  on  spoke  literature. 

As  mentioned  previously,  the  location  within  the  discharge  channel  where  spokes  originate  is  unknown. 
Previous  experiments  have  detected  azimuthal  oscillations  throughout  the  plasma  from  the  anode  out  into  the 
plume.  The  original  work  by  Janes  and  Lowder1  detected  azimuthal  oscillations  with  probes  in  the  discharge 
channel  where  they  even  noted  a  spoke  angle  with  respect  to  the  walls.  However,  significant  differences  exist 
between  their  experimental  apparatus  in  1966  and  the  modern  H6  including  a  longer  discharge  channel  length, 
~  1/2  discharge  voltage  than  nominal  H6  values,  chamber  pressure  two  orders  of  magnitude  higher,  use  of 
filament  cathodes,  radially  inward  magnetic  field  with  over  2  x  larger  peak  value  and  quartz  walls  instead 
of  Boron  Nitride.  More  recent  use  of  probes  downstream  of  the  discharge  channel  exit  plane  by  Chesta23 
and  Sekerak7, 24  observed  spokes  propagating  in  the  plume  at  least  one  channel  diameter  downstream.  Work 
by  McDonald11  with  a  segmented  anode  showed  PSD  frequency  peaks  that  matched  spokes  frequencies 
from  HI  A  PSD.  In  a  time  history  trace  (Figure  9  from  Ref.  11)  local  discharge  current  peaks  are  observed 
propagating  from  segment  to  segment  at  a  velocity  between  1700  and  2500  m/s.  Therefore,  the  presence  of 
azimuthal  oscillations  are  detectable  from  the  anode  out  into  the  plasma  plume  of  HETs  with  theories  and 
experiments  disagreeing  on  where  they  originate. 

In  searching  for  the  origin  for  spokes,  we  will  consider  the  plasma  in  three  different  regions  similar  to 
that  used  in  simulations:25  near-anode  region,  the  ionization/acceleration  region  near  the  channel  exit,  and 
near- field  plasma  plume.  HETs  have  steep  gradients  in  all  relevant  plasma  parameters:  density,  electron 
temperature,  electric  field  and  magnetic  field.  However,  the  largest  gradients  occur  near  the  channel  exit 
where  the  magnetic  field  peaks  and  the  bulk  ionization  and  accelerations  occurs.  The  origin  and  mechanism 
for  spokes  are  unknown,  but  we  will  begin  our  discussion  by  broadly  considering  these  different  mechanisms: 

1.  Plasma  Waves  The  values  for  magnetic  field,  plasma  density,  neutral  density,  plasma  potential  and 
electron  temperature  vary  in  space  throughout  the  discharge  channel.  However,  if  they  are  constant 
in  time,  slowly  varying  on  time-scales  longer  than  spoke  propagation  time-scales  or  have  a  small 
oscillation  amplitude,  then  any  number  of  plasma  waves  can  propagate  that  can  be  described  by  a 
dispersion  relation.  Section  VC  discussed  the  implications  of  this  in  more  detail  as  a  possible  indicator 
of  breathing- mode  damping.  We  can  subdivide  the  plasma  wave  mechanisms  into  three  categories: 

(a)  (Nearly)  Homogeneous  Plasma  Waves.  The  gradients  in  the  near-anode  region  or  near- field 
plasma  plume  may  be  sufficiently  small  such  that  elementary,  homogenous  plasma  waves  may 
exist.  Discussed  in  Section  IVD. 

(b)  Drift  Waves  Driven  by  Gradients  in  the  Plasma.  The  ionization  and  acceleration  region 
near  the  discharge  channel  exit  has  steep  gradients  in  all  plasma  parameters  and  provides  fertile 
ground  for  drift  waves  without  considering  ionization  effects.  A  preliminary  investigation  for  this 
is  conducted  in  Section  IVE  and  a  more  comprehensive  analysis  of  the  plethora  of  drift  dispersion 
relations  (Section  IVE3)  reserved  for  future  work. 

(c)  Azimuthal  Ionization  Instability.  Chesta’s26  and  Escobar’s22  linear  analyses  showed  oscilla¬ 
tions  due  to  ionization  and  Ref.  15  suggests  that  spokes  represent  increased  ion  production  zones. 
These  waves  are  dependent  on  ionization  frequencies  and  neutral  density,  which  are  captured  in 
some  of  the  dispersion  relations  identified  in  Section  IVE3.  Discussion  of  these  mechanisms  are 
reserved  for  future  work. 

2.  Localized  Breathing  Mode  The  same  fundamental  mechanism  may  exist  between  the  spoke  and 
breathing  modes,  where  the  breathing  mode  exhibits  a  uniform  channel  discharge  and  replenishment 
process  while  the  spoke  mode  exhibits  an  azimuthally  local  discharge  and  replenishment  process.  As 
shown  by  several  numerical  simulations27  29  and  in  Fig.  12,  the  breathing  mode  is  a  slow  progression  of 
the  neutral  front  down  the  channel  leading  up  to  a  rapid  ionization  event  that  Barral30  calls  avalanche 
ionization.  During  the  neutral  replenishment  period  within  the  discharge  channel,  slight  perturbations 
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or  azimuthal  non-uniformities  in  electron  density,  electron  temperature  and/or  neutral  density  (all 
related  to  ionization  rate),  cause  a  local  region  in  the  discharge  channel  to  reach  the  avalanche  ioniza¬ 
tion  point  before  adjacent  regions,  thus  causing  a  localized  increase  in  plasma  density.  The  localized 
electron  density  enhancement  travels  azimuthally  in  the  E  x  B  direction  initiating  avalanche  ionization 
in  an  adjacent  region  thus  propagating  the  perturbation.  Extensive  numerical  and  experimental  char¬ 
acterization  of  the  neutral  flow  and  plasma  within  the  discharge  channel  was  performed  on  the  H631 
where  a  primary  ionization  region  of  ^25%  of  the  discharge  channel  length  and  neutral  flow  velocity 
of  250  m/s  can  be  reasonably  inferred.  These  assumptions  yield  a  neutral  refill  time  of  the  ionization 
region  of  ^40  /is,  corresponding  to  a  25-kHz  oscillation,  which  is  in  line  with  the  oscillation  frequencies 
observed.  Breathing  mode  is  typically  lower  frequency  than  spokes,  so  in  breathing  mode  the  neutral 
front  could  be  receding  deeper  into  the  channel.  Detailed  investigation  of  this  mechanism  is  reserved 
for  future  work. 

3.  Wall  Effects.  Recent  research  suggests  that  plasma  contact  with  the  wall  or  potentially  wall  heating 
could  be  related  to  spoke  propagation  and  mode  transition.  Discussed  briefly  and  only  qualitatively  in 
Section  VD. 


B.  H6  Internal  Data 

Internal  measurements  were  made  by  Reid31  on  the  H6  at  300  V  with  20  mg/s  anode  flow  rate  after  less 
than  300  hours  of  operation.  The  magnet  settings  used  were  J/m  =  3. 00 A,  Iqm  =  2.68  A  which  corresponds 
to  Br/B*  =  0.86.  Table  1  shows  the  source  of  internal  data  used  for  calculations  in  the  following  sections. 
These  figures  have  been  reproduced  in  Fig.  14(a)- (e)  of  the  Appendix  for  reference. 


Variable  Source 
rii  Reid31 

Te  Reid31 
Vp  Reid31 

Ez  Reid31 
N  Reid31 

Table  1:  Reference 
in  Fig.  14(a)-(e). 


Figure  7-10,  20  mg/s 

Figure  7-5,  20  mg/s 
Figure  7-18,  20  mg/s 

Figure  7-21,  20  mg/s 
Figure  8-9,  20  mg/s 


Notes 

Blended  solution  from  OML  and  thin  sheath  Langmuir 
probe  measurements 

Langmuir  Probe,  compare  with  Figure  15  from  Hofer10 
Emissive  probe  corrected  with  Te  from  Langmuir  probe, 
compare  with  Figure  15  from  Hofer,10  referenced  to  cathode 
Computed  from  derivative  of  plasma  potential 
Computed  from  1-D  heavy  particle  continuity  analysis 


and  notes  for  plasma  measurements  of  H6  discharge  channel  (internal).  Data  is  shown 


C.  Frequencies 

The  frequencies  of  interest  for  this  investigation  are  5-15  kHz  for  breathing  mode  and  5-30  kHz  for  spokes.  For 
comparison,  the  cyclotron  frequency  is  uici,e  =  \q\  B / Vn^e  and  the  plasma  frequency  is  ujpi^e  =  yjnq2 / (eo m^e), 
where  subscript  i  is  for  ions  and  e  is  for  electrons.  Variables  include  q  for  the  elementary  charge,  B  for  the 
magnitude  of  the  magnetic  field,  m^e  are  the  ion  and  electron  mass,  n  for  the  plasma  density  where  quasi¬ 
neutrality  has  been  assumed  such  that  ni  ~  ne  ~  n,  and  eo  is  the  permittivity  of  free  space.  The  lower 
hybrid  frequency32  is  uJih  —  sj^d^ce  for  ujpi  uJd-  Table  2  shows  these  characteristic  frequencies  at  three 
different  locations  in  the  discharge  channel  z  =  0.25,  1.00  and  1.5  Lchni  based  on  the  data  in  Fig.  14,  but 
can  be  summarized  as:  UJd  <  u  <  uJih  <  upi  <  ujce  <  u)pe . 

D.  Homogeneous  Plasma  Waves 

Here  we  consider  some  simple  drifts  and  homogeneous  plasma  waves  that  propagate  perpendicular  to  mag¬ 
netic  fields  such  as  E  x  B  drift,  electrostatic  ion  cyclotron  waves  and  magnetosonic  waves.  Other  waves 
summarized  by  the  Clemmow-Mullaly- Allis  diagram33  that  can  propagate  perpendicular  to  magnetic  fields 
such  as  ordinary  waves  (O  wave),  extra  ordinary  waves  (X  waves)  and  upper- hybrid  resonance  are  of  higher 
frequency  than  the  10’s  kHz  spoke  modes. 


16 

The  33rd  International  Electric  Propulsion  Conference ,  The  George  Washington  University,  USA 

October  6-10,  2013 

Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


0.25Lchnz  1.00Lchni  1.50  ^chni 


Oscillation 


Ion  cyclotron,  kHz 
Lower  hybrid,  kHz 
Ion  plasma,  MHz 
Electron  cyclotron,  MHz 
Electron  plasma,  GHz 


0.25 

1.7 

1.3 

120 

810 

650 

9.4 

27 

20 

61 

400 

320 

4.6 

13 

9.6 

Table  2:  Representative  frequencies  on  channel  centerline  for  Region  I,  II  and  III  at  0.25,  1.00  and  1.50 
Lchni ,  respectively.  Ion  cyclotron  and  ion  plasma  frequencies  are  for  Xe1+. 


1.  E  x  B  Drift 

There  are  many  single  particle  drifts  to  consider  and  only  a  few  will  be  considered  here;  future  work  will 
investigate  the  VH  and  curvature  drift.  The  E  x  B  drift  velocity  and  6  and  component  are  given  by34 

E  x  B 


z  EzBr  -  ErBz 
ve,b0-  m  +  m 

The  calculated  E  x  B  drift  velocity  distribution  in  the  discharge  channel  is  shown  in  Fig.  14(f)  where  the 
maximum  value  is  over  5  x  106  m/s  at  the  peak  electric  field  and  of  order  105  m/s  within  ±0.2Lchni-  These 
velocities  are  two  to  three  orders  of  magnitude  higher  than  the  typical  spoke  velocity  of  1500  to  2200  m/s. 
It  is  interesting  to  note  the  E  x  B  drift  velocity  is  of  the  same  order  as  the  electron  thermal  velocity  shown 
in  Fig.  14(g). 

2.  Electrostatic  Ion  Cyclotron  Waves 

Similar  to  sound  waves  in  air,  ion  acoustic  waves  follow  a  simple  dispersion  relation 


(16) 

(17) 


uo  =  k\s 


(18) 


where  vs  is  the  ion  acoustic  speed  with  electrons  of  temperature  Te  in  eV 


v 


S  — 


(19) 


As  shown  in  Fig.  14(c)  the  electron  temperature  in  most  of  the  channel  and  the  near  field  plume  is  ^  5  eV, 
except  for  near  the  exit  plane  0.7  <  Lchni  <1.1  where  the  electron  temperature  peaks  at  ~  35  eV.  The  ion 
acoustic  speed  for  Te  =  5  eV  is  2000  m/s  and  for  Te  =  35  eV  is  5000  m/s.  Fig.  14(h)  shows  the  acoustic 
speed,  which  is  close  the  characteristic  speed  of  ^  2000  m/s  for  most  of  the  discharge  channel.  Unfortunately, 
ion  acoustic  waves  propagate  parallel  to  the  magnetic  field  while  spokes  propagate  perpendicularly. 

An  electrostatic  ion  cyclotron  wave33  is  similar  to  an  ion  acoustic  oscillation  except  the  Lorentz  force 
provides  a  restoring  force32  which  yields  a  modification  to  Eq.  18 


UJ 


7  2  2  I  2 

k  v*  + 


Ion  acoustic  waves  can  propagate  nearly  perpendicular  to  B  and  have  a  phase  velocity  of 


(20) 


vp/*  =  \fl2s  +wci/fc2  (21) 

Except  for  the  difference  in  sign  inside  the  radical,  note  the  similarity  to  Eq.  8  with  a  =  2  with  the  ion 
acoustic  speed  as  the  characteristic  velocity  and  the  ion  cyclotron  frequency  as  the  characteristic  frequency. 
Since  spoke  location  is  unknown,  spokes  could  be  related  to  electrostatic  ion  cyclotron  waves  in  the  channel 
near  the  anode  or  in  the  near  field  plume  region  where  Te  ^  5  eV.  A  map  of  the  electrostatic  ion  cyclotron 
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frequencies  is  shown  in  Fig.  14(i)  for  kg  corresponding  to  m  =  5  where  the  frequencies  are  within  the 
expected  range  for  spokes.  Unfortunately,  this  mechanism  only  implies  that  waves  are  perpendicular  to  the 
magnetic  field  (axial,  CW  azimuthal  or  CCW  azimuthal)  and  does  not  force  spokes  to  propagate  in  the 
E  x  B  direction  (CCW  azimuthal). 


3.  Magnetosonic  Waves 


Magnetosonic  waves  are  low  frequency,  electromagnetic  waves  propagating  across  the  magnetic  field.  They 
are  similar  to  acoustic  waves,  but  the  oscillations  are  produced  by  E  x  B  drifts  across  E32  The  phase 
velocity  is 


Vph  = 


(22) 


where  v^  =  B / yj nomin  is  the  Alfven  speed  and  /iq  is  the  permeability  of  free  space.  Fig.  14(j)  shows  the 
magnetosonic  phase  velocity  (which  is  very  close  to  the  Alfven  speed  so  it  is  not  shown  in  Fig.  14),  both  of 
which  are  an  order  of  magnitude  too  large  to  be  azimuthal  spokes. 


E.  Gradient  Drift  Waves 


Drift  waves  are  common  when  spatial  gradients  exist  in  plasma  properties  because  they  provide  free  energy 
from  which  an  instability  can  grow  and  a  wave  can  propagate.  Figs.  14  (a)-(e)  clearly  shows  the  plasma  in 
the  discharge  channel  and  near-field  plume  are  not  uniform  therefore  spatial  gradients  exist,  predominantly 
near  the  channel  exit.  There  are  many  forms  of  drift  waves  that  have  been  studied  from  fusion  research. 
Here  we  start  with  the  most  common  form  of  the  drift  wave  that  arises  from  a  gradient  in  plasma  density  and 
proceed  with  increasing  complexity.  Also  important  is  the  definition  of  gradient  length  scale  for  parameter  g 


Ls/q  — 


(23) 


1.  Density  Gradient  Drift  Waves 

The  most  fundamental  gradient  driven  wave  is  the  density  gradient  which  assumes  no  steady  state  electric 
field  and  constant  magnetic  field,  neither  of  which  are  fundamentally  true  in  a  HET  discharge  channel.  The 
dispersion  relation  for  density  gradient  drift  waves  and  the  diamagnetic  drift  velocity  are33 

lj  =  kgy  Drift 

_Te  1  dne  _  Te 
VDHfi  ~  Bne  dz  ~ 

Drift  waves  are  caused  by  density  perturbations  from  ion  E  x  B  drift  where  ions  ”  slosh”  back  and  forth  in 
the  gradient  (axial)  direction,  with  the  phase  velocity  in  the  Vne  x  B  direction.34  Therefore,  no  particles 
or  energy  are  moving  in  the  theta  direction.  In  the  H6,  Vne  is  in  the  +z  direction  from  the  anode  to 
L  <  0.8 Lchni  and  is  in  the  —  z  direction  for  L  >  Lchni  from  Fig.  14(a).  Drift  waves  would  be  in  the  E  x  B 
for  L  <  O.SLchni  an(4  in  the  —E  x  B  for  L  >  Lchni •  Spokes  always  propagate  in  the  E  x  B  direction,  so  we 
focus  on  the  density  gradient  upstream  from  the  ionization  zone  for  density  gradient  waves  to  be  the  source 
of  spoke  motion. 

The  measured  temperate  in  the  discharge  channel  varies  from  5  to  35  eV  as  shown  in  Fig.  14(c).  Although 
the  magnetic  field  magnitude  and  profile  cannot  be  shown,  it  is  of  order  100’s  G  similar  to  the  SPT-100. 
Therefore,  the  Te/B  term  in  Eq.  25  is  of  order  103  which  is  the  same  as  the  spoke  velocity.  The  density 
length  scale  term  L\7ne  must  be  of  order  unity  for  v Drift  to  be  in  the  spoke  velocity  range  of  1500-2200  m/s. 
The  experimental  data  from  Reid  in  Fig.  14(a)  for  ion  density  on  channel  centerline  shows  an  unexpected 
multi-peaked  feature  that  is  likely  experimental  error  so  the  axial  derivative  of  density  is  smoothed  with  a 
0.0 6Lchni  moving  average  filter.  The  density  gradient  length  scale  is  of  order  Lvne  ^  10-2  causing  the  drift 
velocity  to  be  approximately  two  orders  of  magnitude  too  large.  Discounting  the  unusual  multi-peaked  ion 
density  yields  a  similar  result  where  the  density  increases  from  1  x  1018  to  4  x  1018#/m3/  in  0.2 Lchni  which 
corresponds  to  dn/dz  =  3.6  x  1020#/m3/m  and  v Drift  ~  105.  The  drift  velocity  is  shown  in  Fig.  14(k) 
for  the  entire  discharge  channel  and  is  larger  than  the  spoke  or  characteristic  velocities  we  are  seeking.  It  is 
interesting  to  note  that  this  velocity  is  of  the  same  order  or  lower  than  E  x  B  drift  velocity. 
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(24) 

(25) 


2.  Collisional  Density  Gradient  Drift  Waves 


The  next  step  is  to  include  collisional  drift  waves,  which  in  the  limit  of  small  propagation  along  the  B-field 
lines  for  a  simple  slab  (cartesian  coordinates)  is34 

_  kevprift 

^  _  1  +  kgVg/uJci 

The  phase  velocity  for  Eq.  26  is  shown  in  Fig.  14(1)  for  kg  corresponding  to  m  =  5,  which  yields  results 
lower  than  the  expected  spoke  velocity  by  as  little  as  50%  to  as  much  as  an  order  of  magnitude.  Noting 
that  kgvl/ujci  ^  1?  Eq.  26  reduces  to  u  ^  ujcij  (kgL^/nf)^  which  is  not  the  same  functional  form  as  Eq. 
6.  Therefore,  the  spokes  are  unlikely  caused  by  density  gradient  drift  waves  or  collisional  drift  waves,  in  a 
simplified  slab  geometry. 

3.  Future  Investigation  for  Drift  Waves 

The  body  of  literature  on  dispersion  relations  for  various  drift  waves  is  voluminous  and  many  other  appropri¬ 
ate  relations  should  be  considered  in  the  same  manner  as  above.  Future  work  will  investigate  the  collisional 
drift  instability  in  cylindrical  coordinates  as  developed  by  Ellis35  which  was  successfully  used  by  Jorns36  to 
identify  drift  waves  related  to  cathode  oscillations  in  magnetically  shielded  thrusters.  Additionally,  a  drift 
relation  shown  by  Esipchuk37  that  accounts  for  density  and  magnetic  field  gradients  has  been  used  to  repro¬ 
duce  azimuthal  oscillations  in  the  10’s  kHz  in  the  near-field  plume  of  an  SPT-100.21  Kapulkin38  developed 
a  similar  dispersion  relation  analysis  for  the  near-anode  region  of  an  SPT-100,  which  is  also  very  applicable 
to  the  H6  conditions  tested  here.  Frias  recently  developed  two  new  dispersion  relations;  one  accounts  for 
electron  flow  compressibility  and  the  other  includes  temperature  oscillations  and  gradients.39,40  Finally,  the 
linearized  2-D  axial- azimuthal  models  by  Chesta26  and  Escobar22  accounts  for  ionization  and  neutral  density 
should  be  investigated.  Applying  all  of  these  dispersion  relations  to  the  plasma  parameters  described  above 
and  comparing  to  the  empirical  dispersion  analysis  of  Section  III  could  yield  valuable  into  the  origin  and 
nature  of  azimuthal  spokes. 


(26) 


V.  Discussion 


A.  Spoke  Velocities 

The  manual  method,  correlation  method  and  dispersion  method  all  yield  very  similar  results  for  the  spoke 
velocity  building  confidence  in  all  three  techniques.  In  the  future,  spoke  velocities  will  be  reported  routinely 
as  an  additional  metric  with  high-speed  imaging  data. 

The  inverse  relation  of  spoke  velocity  to  magnetic  field  is  an  interesting  trend  that  should  help  identify 
the  spoke  mechanism  even  though  the  spoke  velocity  is  only  weakly  dependent  on  magnetic  field  ( vsp  oc  B~ & 
where  0.25  <  (3  <  0.5).  A  previous  study6  showed  similar  results,  but  the  spoke  velocity  for  each  spoke 
order  m  was  calculated  from  ui^/kg  (which  effectively  gives  the  characteristic  velocity)  instead  of  the  more 
rigorous  techniques  used  here.  The  results  in  Fig.  7(a)  can  be  compared  with  Figure  5  from  Ref.  6,  which 
shows  the  same  trend  of  decreasing  spoke  velocity  up  to  point.  However,  the  same  H6  thruster  was  used  and 
lower  spoke  orders  were  seen  to  dominate  in  that  study  for  unexplained  reasons.  In  contrast,  the  early  work 
of  Lomas41  found  that  spoke  velocity  increased  with  magnetic  field,  which  is  contrary  to  what  was  observed 
by  McDonald6  and  in  this  investigation. 

The  E  x  B  drift  velocity  also  appears  to  have  a  simple  1/B  dependence  at  first  glance.  However,  noting 
that  E  is  related  to  B  from  E  =  77  (l  +  fl2)  j  where  g  is  the  plasma  resistivity,  j  is  the  discharge  current 
density  and  De  is  the  Hall  parameter,  means  the  dependence  of  vexB  cm  B  is  actually  more  complicated. 
Spokes  propagate  in  the  E  x  B  direction,  so  they  may  be  related  to  electron  or  possibly  even  ion  E  x  B  drift. 
The  Larmor  radius  for  electrons  is  less  than  0.05 Lchni  using  their  thermal  velocity  so  they  are  magnetized 
as  expected.  The  ion  Larmor  radius  for  +1  ions  using  the  ion  velocity  calculated  from  plasma  potential 
and  energy  conservation  is  greater  than  10 Lchni  for  most  of  the  discharge  channel  except  for  the  ionization 
region  before  they  are  accelerated  by  the  large  electric  field.  Note  the  Larmor  radius  for  ions  decreases  for 
higher  charge  states,  which  are  known  to  exist.  If  the  ionization  region  is  sufficiently  offset  upstream  from 
the  acceleration  region,  then  a  region  of  large  ion  density  with  relatively  low  ion  velocity  may  exist  ( <  5  x  103 
m/s,  of  same  order  as  spoke  velocity).  In  this  region,  the  ions  may  be  able  to  undergo  some  azimuthal  motion 
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before  being  accelerated  downstream,  although  unlikely  completing  a  full  cyclotron  orbit.  This  motion  could 
contribute  to  azimuthal  spokes. 

The  observation  that  ion  acoustic  speed  matches 
the  characteristic  speed  from  the  dispersion  method 
is  encouraging  because  vs  commonly  appears  in 
waves  such  as  the  electrostatic  ion  cyclotron  wave 
and  arises  prominently  in  drift  waves.  For  Br/B*  = 

0.86,  the  spoke  velocity  is  1543  m/s  from  Fig.  11  and 
the  characteristic  velocities  are  2188  and  1848  m/s 
for  a  =  1  and  2  from  Fig.  10(a).  Fig.  12  shows  a 
comparison  of  those  velocities  with  the  channel  cen¬ 
terline  ion  acoustic  velocity  from  Fig.  14(h)  and  the 
critical  ionization  velocity.  Janes  and  Lowder  sug¬ 
gested  that  the  spokes  may  be  related  to  the  criti¬ 
cal  ionization  velocity1  first  proposed  by  Alfven  and 
Ref.  6  contains  a  good  discussion  of  the  phenomenon 
and  the  implication  to  HETs.  The  critical  ionization 
velocity  for  Xenon  is  4200  m/s  and  shown  in  Fig.  12, 
which  is  the  same  order  of  magnitude  as  the  spoke 
velocity,  but  is  still  over  twice  the  value.  The  ion 
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Figure  12:  Ion  acoustic  speed  on  channel  center  line  for 
Br/B*  —  0.86  from  Fig.  14(h)  smoothed  by  a  0.0375 
z/Lchni  moving  average  filter.  The  critical  ionization 
velocity,  spoke  velocity  and  characteristic  velocities  for 
a  =  1,  2  are  shown  for  comparison. 


acoustic  speed  matches  the  characteristic  speed  for  a  =  2  better  than  a  =  1,  particularly  for  z/Lchni  <  0.7. 
The  spoke  velocity  is  lower  than  the  ion  acoustic  speed  in  the  near-anode  region,  but  is  similar  for  the  near 
field  plume  region.  The  similarity  between  ion  acoustic  speed  and  spoke  velocity  or  characteristic  velocity  is 
unlikely  a  coincidence  and  will  be  further  investigated  with  the  more  detailed  dispersion  relations  discussed 
in  Section  IVE3. 


B.  High  B-field  Mode 

The  discussion  in  Section  II  identified  two  different  modes  for  HET  operation:  global  oscillation  mode  and 
local  oscillation  mode.  However,  the  plasma  oscillations  in  the  discharge  channel  show  different  characteristics 
for  high  magnetic  fields  that  may  be  indicative  of  yet  another  different  operational  mode.  Fig.  2  shows  that 
as  the  magnetic  field  is  increased,  the  spoke  order  increases,  the  spoke  duration  decreases  and  low-amplitude, 
entire  channel  oscillations  occur.  Additionally,  Figures  12  and  13  of  Ref.  7  show  the  discharge  current  RMS 
increases  at  high  magnetic  field  settings  for  all  conditions  except  300  V,  25.2  mg/s.  Finally,  Fig.  11  shows 
that  spoke  velocity  becomes  independent  of  magnetic  field  for  Br/B*  >  1.  The  transition  to  this  mode 

is  not  as  sharp  as  the  distinct  transition  between  global  and  local  mode.  A  qualitative  explanation  can  be 
offered  if  collisional  drift  waves  are  responsible  for  azimuthal  spokes.  In  discussing  collisional  drift  waves, 
Bellan  notes34  ” plasmas  with  strong  magnetic  fields  tend  to  have  turbulent,  short  perpendicular  wavelength 
drift  waves,  whereas  plasmas  with  weak  magnetic  fields  have  coherent,  long  perpendicular  wavelength  drift 
waves.”  This  could  indicate  that  the  high  magnetic  field  settings  are  a  turbulence  dominated  oscillation 
mode. 


C.  Breathing  Mode  Damping 

The  investigation  on  HET  mode  transitions7  showed  that  spokes  were  present  and  dominated  plasma  oscilla¬ 
tions  when  the  global  oscillation  ceased.  We  can  postulate  that  the  global  oscillation  mode  is  the  well-studied 
breathing  mode  and  the  transition  to  local  oscillation  mode  represents  a  damping  or  cessation  of  the  breath¬ 
ing  mode  allowing  the  azimuthal  spokes  to  propagate.  The  breathing  mode  has  been  numerically  modeled  by 
Fife,27  Boeuf,28  Barral29  and  Hara,42  however  Choueri21  pointed  out  in  2001  that  stability  criteria  has  yet  to 
be  theoretically  developed.  Recently,  Barral43  published  the  first  known  stability  analysis  for  the  breathing 
mode  using  linear  and  non-linear  techniques,  but  further  validation  is  necessary  to  confirm  the  criteria. 

All  of  the  preceding  analysis  on  spoke  propagation  and  dispersion  relations  assumes  the  plasma  properties 
are  constant  in  time  within  the  discharge  channel  allowing  spokes  to  propagate.  In  the  H6,  a  typical  spoke 
will  require  ^  300 /is  to  travel  around  the  entire  discharge  channel.  The  1-D  fluid  simulation  of  an  SPT-100 
developed  by  Barral29  shows  that  during  a  ~  100/is  breathing  mode  cycle  the  plasma  properties  can  change 
significantly.  Fig.  13  for  the  plasma  properties  on  discharge  channel  centerline  shows  that  neutral  density  can 
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change  by  a  factor  of  2-3,  plasma  density  can  change  up  to  an  order  of  magnitude,  and  electron  temperature 
can  change  by  a  factor  of  2  (from  Figure  5  of  Ref.  29,  not  reproduced  here).  The  simulation  calculated  a 
slow  progression  of  the  neutral  front  down  the  discharge  channel  building  up  to  avalanche  ionization.  Fig. 
13  shows  that  once  avalanche  ionization  occurs  at  ~  50  /is  and  again  at  ~  150  /is,  the  plasma  density  peaks 
and  the  neutral  front  recedes  with  a  velocity  in  the  range  of  0.5  —  5  x  103  m/s.  Spokes  may  not  be  able  to 
propagate  since  the  breathing  mode  period  is  of  the  same  order  or  smaller  than  the  time  required  for  spoke 
propagation  and  the  recession  rate  can  be  greater  than  the  spoke  velocity  of  ^  1700  m/s. 

A  related  paper  at  this  conference  by  Hara44  uses  a  hybrid 
direct-kinetic  simulation  to  vary  magnetic  field  and  induce  a 
mode  transition  in  an  SPT-100  type  thruster.  There  it  is  found 
that  the  breathing  mode  is  damped  in  a  region  of  optimal  per¬ 
formance  similar  to  the  experimental  results  of  Ref.  7.  A  future 
publication  will  investigate  the  breathing  mode  frequency  vari¬ 
ation  and  transition  criteria  using  a  fluid  model.  The  results 
from  Ref.  7  will  be  combined  with  the  1-D  transient  fluid  code 
in  Ref.  29  to  experimentally  and  computationally  investigate 
breathing  mode  stability  and  assess  whether  HET  mode  tran¬ 
sitions  are  indeed  breathing  mode  damping. 


D.  Wall  Considerations 

A  HET  does  not  naturally  have  a  radial  restoring  force  like 
a  magnetron  so  the  electrons  likely  interact  with  the  wall  in 
some  way  (magnetic  mirror  or  radial  electric  field)  to  propa¬ 
gate  their  motion  to  make  HETs  closed-drift  devices.  Recent 
research  shows  that  plasma  contact  with  the  wall  or  potential 
wall  heating  could  be  related  to  spoke  propagation  and  mode 
transition.  The  transition  between  modes  was  noted  to  be  af¬ 
fected  by  time  at  the  given  operating  condition,  which  could 
be  related  to  wall  heating  and  thermal  equilibrium,7  although 
wall  temperatures  were  not  measured  in  that  investigation.  Ad¬ 
ditionally,  while  investigating  magnetically  shielded  thrusters 
where  the  plasma  is  not  in  contact  with  the  walls,  Jorns36  did 
not  observe  spokes.  The  magnetically  shielded  H6  showed  sim¬ 
ilar  performance  to  the  unshielded  H610  suggesting  the  spokes 
are  only  prominent  in  non- magnetically  shielded  thrusters.  Additionally,  thrusters  with  partial  ceramic  walls 
(boron  nitride  near  the  channel  exit /ionization  zone  and  conducting  walls  near  the  anode)  may  show  different 
spoke  behaviors  if  wall  effects  are  important.  Evidence  of  these  different  behaviors  will  be  investigated  in 
future  work. 

While  investigating  different  wall  materials  in  the  SPT-100,  Gascon45  observed  mode  transitions  similar 
to  those  observed  in  Ref.  7.  Gascon  did  not  have  the  ability  to  observe  spokes  in  that  investigation 
so  we  cannot  confirm  that  spokes  were  present  in  the  low  discharge  current  mode.  The  accompanying 
theory  provided  by  Barral46  suggested  Space  Charge  Saturation,  when  the  wall  secondary  electron  emission 
coefficient  approaches  unity,  was  responsible  for  the  mode  transition.  If  this  is  correct,  it  implies  a  wall 
related  phenomenon  for  mode  transition  and  hence  the  onset  or  disappearance  of  spokes. 

Even  though  we  do  not  have  a  proven  theory  for  spoke  formation  and  propagation,  we  can  speculate 
that  the  electron  E  x  B  drift  motion  inherent  in  HET  discharge  channels  is  related  to  spokes  since  spokes 
always  propagate  in  the  E  x  B  direction.  As  was  discussed  by  King47  based  on  single  particle  motion,  a 
radial  electric  field  at  the  outer  wall  is  required  to  keep  the  electrons  on  a  circular  path  in  the  discharge 
channel  due  to  the  cylindrical  geometry.  Katzb  pointed  out  that  the  sinusoidal  motion  suggested  by  King 
is  likely  not  accurate  because  the  electrons  would  spend  more  time  at  the  extrema  (inner  and  outer  walls) 
causing  the  electron  density  to  peak  off-channel  centerline  which  is  not  seen  in  Fig.  14(a).  We  can  calculate 
the  electric  field  required  to  keep  the  electrons  in  a  circular  path  following  the  discharge  channel  walls  from 
simple  centripetal  acceleration.  Considering  the  E  x  B  drift  velocity  and  the  electron  thermal  velocity,  the 
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Figure  13:  Neutral  density  (top)  and 
plasma  density  (bottom)  from  1-D  fluid 
simulation  of  an  SPT-100  on  channel  cen¬ 
terline.  Reproduced  from  Figure  3(b)  and 
4(b)  of  Ref.  29. 
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range  of  azimuthal  electron  velocities  v#  to  consider  is  1  x  104  to  5  x  106  m/s.  Using  the  electron  equation  of 
motion,  the  radial  electric  field  Er  =  (me / q)^ / Router)  required  for  the  electron  to  maintain  circular  motion 
at  the  outer  wall  is  6  /iV/mm  to  1.5  V/mm.  This  is  a  small  electric  field  that  could  be  established  by  sheaths 
at  the  walls  or  internal  to  the  plasma  in  order  to  maintain  quasineutrality.  Therefore,  the  electrons  are  likely 
in  nearly  circular  motion  around  the  discharge  channel.  A  radial  electric  field  at  the  walls  is  required  to 
reflect  electrons  back  into  the  channel  per  King  or  is  required  to  keep  electrons  in  a  circular  motion  around 
the  discharge  channel.  If  this  field  cannot  be  established  or  becomes  perturbed,  then  the  electron  motion 
around  the  discharge  channel  could  be  disrupted  which  could  disrupt  spoke  propagation  assuming  they  are 
related. 


VI.  Conclusions 

The  spoke  velocity  has  been  determined  using  three  methods  with  similar  results:  manual  fitting  of 
diagonal  lines  on  the  spoke  surface,  linear  cross-correlation  between  azimuthal  locations  and  an  approximated 
dispersion  relation.  The  spoke  velocity  for  three  discharge  voltages  (300,  400  and  450  V)  and  three  anode 
mass  flow  rates  (14.7,  19.5  and  25.2  mg/s)  yielded  spoke  velocities  between  1500  and  2200  m/s  across  a  range 
of  normalized  magnetic  field  settings.  The  spoke  velocity  was  inversely  dependent  on  magnetic  field  strength 
for  Br/B*  <  1  and  asymptoted  to  1600-1700  m/s  for  Br/B *  >  1  for  all  conditions  except  300  V,  14.7  mg/s. 
The  dispersion  plots  for  probes  1.5  Rchni  downstream  displayed  similarities  to  the  dispersion  plots  calculated 
from  HIA  reinforcing  previous  results  that  oscillations  observed  with  imaging  of  the  discharge  channel  are 
correlated  to  oscillations  in  the  plasma  plume.  Spoke  velocity  from  the  probe  delay  method  was  ~  30% 
higher  for  300  V,  19.5  mg/s  and  for  400  V,  19.5  mg/s  showed  an  unexplained  increase  for  Br/B *  >  0.9. 

The  mechanism  for  spoke  formation  is  unknown  although  we  consider  homogeneous  plasma  waves,  gra¬ 
dient  driven  drifts,  ionization  related  instabilities  or  wall-effects.  It  is  unknown  whether  spokes  originate 
in  the  near-anode,  ionization/acceleration  region  or  the  near- field  plume  so  the  plasma  properties  from 
previous  internal  measurements  at  Br/B*  =  0.86  are  used  to  investigate  various  dispersion  relations. 
The  10’s  kHz  oscillations  associated  with  spokes  are  compared  to  standard  plasma  oscillations  showing 
uoci  <  u  <  u)ih  <  upi  <  Cc ;ce  <  tJpe-  The  velocities  and  frequencies  are  compared  to  standard  plasma  waves 
such  as  E  x  B  drift,  electrostatic  ion  cyclotron,  magnetosonic  and  various  drift  waves.  The  E  x  B  drift 
is  two  to  three  orders  of  magnitude  larger  than  the  spoke  velocity  and  is  of  the  same  magnitude  as  the 
electron  thermal  velocity  in  agreement  with  previous  results.  Electrostatic  ion  cyclotron  waves  are  of  the 
same  frequency  as  spokes.  The  empirically  approximated  dispersion  relation  of  u)a  =  —  u“h  where 

a  >  1  yielded  a  characteristic  velocity  that  matched  the  ion  acoustic  speed  for  ~  5  eV  electrons  that  exist  in 
the  near-anode  and  near-field  plume  regions  of  the  discharge  channel.  Simple,  cartesian  coordinate  density 
gradient  drift  waves  have  a  diamagnetic  drift  velocity  that  is  ^2  orders  of  magnitude  higher  and  simple  slab 
collisional  drift  waves  phase  velocities  are  too  low.  While  a  good  start,  these  simple  dispersion  relations  do 
not  cover  the  extensive  body  of  literature  on  drift  waves,  particularly  those  that  have  been  derived  for  HETs 
recently,  and  a  more  extensive  comparison  effort  is  underway  that  will  be  presented  in  future  publications. 

Recently  characterized  plasma  oscillations  in  mode  transitions  have  been  shown  to  greatly  affect  thruster 
performance.  Thrust-to-power  is  maximized  in  the  mode  where  azimuthal  spokes  are  present  and  begins 
to  decrease  by  up  to  25%  when  absent.  Although  spokes  are  not  identified  as  the  cause  of  the  increased 
performance,  their  association  suggests  that  investigating  the  underlying  mechanism  of  spokes  will  benefit 
thruster  operation. 
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Appendix 


(a)  Ion  Density 
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(k)  Density  Gradient  Drift  Velocity  m/s  xlO5 
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Figure  14:  Internal  measurements  of  the  H6  discharge  channel  from  Reid31  as  discussed  in  Table  1  for  (a) 
ion  density,  (b)  neutral  density,  (c)  electron  temperature  (d)  plasma  potential,  and  (e)  axial  electric  field. 
Calculated  velocities  include  (f)  E  x  B  drift  velocity,  (g)  electron  thermal  velocity,  (h)  ion  acoustic  velocity 
and  (j)  magnetosonic  velocity,  (k)  density  gradient  drift  velocity  and  (1)  and  collisional  drift  velocity,  (i)  is 
electrostatic  ion  cyclotron  frequency. 
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Introduction 


Why  are  we  doing  this  work? 


—  Continued  examination  of  alternative  Hall  effect 
thruster  propellants:  Krypton 

—  Interest  in  effects  of  test  facility  environments  on  Hall 
effect  thruster  operation 

—  Fundamental  differences  attributable  to  propellant 
selection:  Kr  vs  Xe 

Why  use  laser  induced  Kr  II  fluorescence  (LIF)? 


—  Non-intrusive  measurements  via  LIF  producing  ion 
acceleration  profiles  and  velocity  distributions 


—  Low  mass  of  Kr  will  increase  diffusion  from  facility 
into  discharge  (~25%  increase;  84  amu  vs  131  amu  Xe) 

What  information  do  we  hope  to  extract? 

—  Continued  optimization  of  Kr  in  a  Hall  effect  thruster 

—  Changes  in  thruster  operation  due  to  Kr  ingestion 
relative  to  Xe 
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Comparison  of  Xenon  to  Krypton  Properties 


Property 

Units 

Xe 

Kr 

Atomic  Mass 

amu 

131.3 

83.8 

1st  Ionization  Energy 

eV 

12.1 

14.0 

2nd  Ionization  Energy 

eV 

21 

24 

3rd  Jonization  Energy 

eV 

32 

37 

Atmospheric  Concentration 

ppb 

87 

1000 

Stable  Isotopes 

9 

6 

Odd  Isotopes 

2 

1 

Critical  Pressure 

MPa 

5.84 

5.50 

Critical  Temperature 

K 

290 

209 

Boiling  Point  (1  atm) 

K 

161 

120 

Krypton  Similar  to  Xenon 

—  Noble  gases,  not  widely  separated 

—  Physical  properties  are  similar 

—  Electrons  bound  more  tightly 

•  Less  electron  shielding 

•  Higher  ionization  energies,  +15% 

—  More  ideal  gas  behavior 

•  Less  compressible 

•  More  difficult  pressurized  storage 


Krypton  vs  Condensibles 

—  Kr  is  nearly  Xe  “drop  in”  replacement 

•  No  changes  required  to  thruster 

•  Minimal  changes  required  to  propellant 
management 

•  No  increased  likelihood  of  S/C 
contamination 

—  Compare  to  Bismuth... 

•  Thruster  design  modifications 

•  High  temperature  propellant  distribution 

•  Significant  likelihood  of  S/C 
contamination 


Busek  BHT-600 


IEPC-2009-115,  Nakles  etal. 
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Power  Hall  Effect  Thruster 
and  Test  Facility 


AFRL  Chamber  6,  Edwards  AFB 

•  Cryogenic  pumping 

•  5  independent  pumps 

•  32,000  I/s 

•  Stainless  Steel  chamber  walls 

•  1.8  m  diameter,  3.0  m  length 

•  Graphite  beam  dump 


Minimum  Pressure:  1.2x1  O'5  Torr 
Maximum  Pressure:  3.2x1  O'5  Torr 

Both  pressures  are  below 
accepted  5x1  O'5  Torr  upper 
limited  for  long  duration  HET 
testing 


[BUSEtO 


Medium  Power  Laboratory  HET 

•  Busek  BHT-600  operating  at  600  W 

•  Single  operating  condition  examined 

•  Previously  measured  xenon  acceleration 

•  Performance  characterized 
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Comparison  of  Global  Parameters 


Krypton  Operating  Parameters 


Kr  Anode  Flow 

25.5  seem 
(1.59  mg/s) 

Kr  Cathode  Flow 

1.5  seem 
(94  /ig/s) 

Anode  Potential 

300  V 

Anode  Current 

1.73  A 

Inner  Coil  Current 

1.75  A 

Outer  Coil  Current 

1.75  A 

Keeper  Current 

0.5  A 

Heater  Current 

3.0  A 

Thrust 

22.4  mN 

Anode  Efficiency 

31% 

Specific  Impulse 

1440  s 

Thrust  measurements  from  Nakles, 
Larson,  and  Harg us,  IEPC-2009-118 


Xenon  Operating  Parameters 


Xe  Anode  Flow 

25.5  seem 
(2.5  mg/s) 

Xe  Cathode  Flow 

1.5  seem 
(147  fig/s) 

Anode  Potential 

300  V 

Anode  Current 

1.93  A 

Inner  Coil  Current 

1.75  A 

Outer  Coil  Current 

1.75  A 

Keeper  Current 

0.5  A 

Heater  Current 

3.0  A 

Thrust 

35.8  mN 

Anode  Efficiency 

44% 

Specific  Impulse 

1460  s 

Xenon  LIF  measurements  from 
Nakles  and  Harg  us,  JPP,  vol.  27, 
no.  4,  2011. 
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Laser-Induced  Fluorescence  Apparatus 
and  Regions  of  Detailed  Examination 


Cfurrt* 

Ml 


More 


Thruftftf 


Single  axis  LIF  with  interior  optical  axis  into  thruster 

—  Similar  to  previous  LIF  apparatus  (Hargus  &  Nakles,  IEEE  Plas  Trans  2008) 

Single  region  of  detailed  examination 

—  Centerline  of  acceleration  channel 


—  Good  measurement  of... 


—  Energy  deposition  (acceleration,  implied  electric  fields) 

—  Ionization  processes  (regions  of  ion  creation) 


—  Compare  to  historical  xenon  data 

—  Same  thruster! 

—  Same  applied  acceleration  potential 
—  Same  magnetic  field 

—  Result:  Very  different  plasma. . . 


Thruster 


Goal:  Compare  Kr  and  Xe 
acceleration  in  the  same  HET 
at  similar  conditions  and 
identify  differences  due  to 
elevated  background  pressure 


Acceleration 
Channel 


Body 


Acceleration 
Channel 


Sample 
Locations 

/ 

Coordinate 
System 
'Ongat  (0.0.0) 
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General  Trends  of  Krypton 
Operation  (Low  Facility  Pressure) 


_ 


*  Xenon  and  krypton  operation  initially  show  very  different  electric  fields 

—  Xenon  (black)  shows  a  higher  field,  peaking  within  thruster  channel 
—  Krypton  (green)  exhibits  a  lower,  broader  E  field,  peaking  at  the  exit  plane 

*  Response  to  thruster  input  variables  (e.g.  B  field,  applied  discharge  potential)  is  similar 

—  Raising  B  field  pushes  peak  field  toward  anode 

—  Electric  field  profile  of  krypton  approaches  that  of  xenon  at  increased  potential  and  B  field 
—  Increasing  discharge  potential  increases  the  electric  field 

*  Electric  field  profile  of  krypton  approaches  that  of  xenon  at  increased  potential  and  B  field 

—  While  krypton  operation  is  not  optimized,  it  appears  similar  to  that  of  xenon  with  respect  to  ion  acceleration! 
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Krypton  Most  Probable  Ion  Velocities  at 
Low  and  High  Facility  Pressures 
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Most  Probable  Velocity  (MPV)  of  the  ions 

—  Peak  of  velocity  distribution 

—  Simple  metric  with  low  ambiguity,  readily  identifiable  even  with  noisy  data 

—  Used  to  produce  E  fields  in  previous  chart 

Low  and  High  facility  pressure  operation 

—  Low  pressure  case  of  minimum  facility  pressure  1 .2E-5  Torr  (Kr  corrected) 

—  High  pressure  case  with  40  seem  Kr  injection  at  3.2E-5  Torr  (Kr  corrected) 

—  Previous  study  showed  dramatic  difference  in  internal  ion  acceleration  (Ref.  1) 

Krypton  300  V  case  (low  vs  high  pressure) 

—  Virtually  no  difference  in  MPV 

—  Velocity  distributions  did  broaden  some,  but  peaks  unchanged  (next  page) 

Krypton  350  V  case  (low  vs  high  pressure) 

—  MPV  increased  slightly,  as  much  as  1  km/s 

—  Velocity  distributions  also  broadened,  but  less  so  than  300  V  case 

Unexpectedly  krypton  operation  does  not  produce  the  substantial 
changes  shown  by  xenon  for  the  same  conditions 
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Axial  Velocity  (m/s) 
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Ion  Population  Fraction  Ion  Population  Fraction 


Kr  Integrated  Populations  within 

Thruster  Channel 


0.1 


0.1 


0.01 


0.01 


[on  Energy  (eV) 


Ion  Energy  (eV) 


Integrated  ion  energy  populations 

—  Converted  VDFs  to  energy  distributions  (1/2  mv2) 

—  Running  integration  between  0  and  1  (100%) 

Ion  acceleration  and  late  term  ionization 

—  Late  term  ionization  prevalent  due  to  low  propellant  utilization 

—  Ionization  creates  low  energy  ion  populations 

—  Monotonic  spacing  of  lines  expected  if  ionization  rate  remains 
constant 

—  Crossing  of  positional  integrations  indicates  substantial 
increases  of  ionization  near  exit  plane  (both  cases) 

Differences  between  low  and  high  pressure  cases 

—  Increase  of  50%  low  energy  ions,  apparently  locally  ionized 

—  Likely  evidence  of  ingestion 

—  Subtle  differences 

—  Only  at  population  fractions  less  than  0.2  diverge 

—  Similar  general  behavior 
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Ion  Population  Fraction  |on  Popu|aiton  Fraction 


Kr  Integrated  Populations  within 
Thruster  Channel  (Ref  #1) 


Ion  Energy  (eV) 


Integrated  ion  energy  populations  for  xenon 

—  Taken  from  Ref  #  1 

—  Same  flow  and  thruster  inputs,  similar  back  pressures 

—  3x  lower  integration  time  degrades  SNR  ~2x 

—  Data  below  -0.03  population  faction  questionable 

—  Xe  data  substantially  noisier  near  exit  plane 

—  Dashed  lines  for  LARGE  uncertainty  below  0.1 
population  fraction 

—  Especially  in  high  pressure  case 

Effects  of  elevated  pressure  for  xenon 

—  Xe  cases  appear  to  only  show  overlapping  integration 
curves  for  high  pressure  case 

—  Limited  basis  of  comparison  due  to  low  SNR 

Comparison  to  krypton  case 

—  Xenon  shows  less  population  overlap,  less  ingestion 

—  Xenon  shows  overlap  for  only  high  pressure  case 

—  SNR  precludes  definitive  connection,  but  provides 
corroborating  evidence  that  ingestion  is  affecting  the 
ion  population  and  energy  distribution 
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Density  (m 


4  J  Magnitudes  of  Exit  Plane  Densities 

for  BHT-600  Thruster  (Kr  &  Xe) 


What  are  the  relative  densities  of  the 
neutrals  and  ions  exiting  the  thruster 
and  the  background  facility  density? 


Exiting  neutral  density  calculated  from  Xe 
neutral  velocity  measurements,  approximated 
by  2x  room  temperature  sonic  velocity,  and 
propellant  utilization  fraction 


Exiting  ion  density  calculated  from  Xe  &  Kr 
ion  velocity ,  and  propellant  utilization  fraction 


Facility  density  calculated  as  ideal  gas 
at  300  K.  Note,  cryogenic  systems 
may  reduce  temperature  and  facility 
densities  may  be  substantially  higher 


0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

Propellant  Utilization 

*  Krypton  cases  -  facility  pressure  has  no  substantial  effect 

—  Exit  neutral  density  » facility  neutral  density  (both  cases) 

*  Xenon  cases  -  facility  pressure  has  substantial  effect 

—  High  pressure  case:  exit  neutral  density  ~  facility  neutral  density 
—  Low  pressure  case:  exit  neutral  density  » facility  neutral  density 


Performance  measurements  of  krypton 
imply  propellant  utilization  of  0.8x  that  of 
xenon  (taken  to  be  0.9),  yielding  krypton 
propellant  utilization  of  ~0. 12 


Relative  density  of  neutrals  at  exit 
plane  (thruster  vs  environment) 
appears  to  produce  substantial 
changes  in  discharge  dynamics 
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Summary  and  Conclusions 
and  Future  Work 


*  Krypton  thruster  operation  responds  to  control  variables  similarly 
to  xenon 

—  Raising  the  magnetic  field,  moves  acceleration  toward  anode 

—  Increased  facility  pressure  moves  acceleration  toward  anode 

*  Krypton  shows  substantially  less  susceptibility  to  facility  neutral 
pressure 

—  Xenon  interior  acceleration  changes  with  facility  pressure  dramatically 

—  Krypton  interior  acceleration  shows  only  small  changes  in  shape  of  velocity 
distribution 

*  Neutral  densities  exiting  thruster  and  of  facility  appear  to  be  critical 

—  Krypton  operation  has  substantially  lower  propellant  utilization  (0.7),  hence 
density  of  neutrals  exiting  thruster  remain  far  above  facility  background 

—  Xenon  operation  has  good  propellant  utilization  (0.9),  hence  density  of  neutrals 
exiting  thruster  approach  that  of  high  facility  pressure  case 

—  Low  facility  pressure  case  has  neutral  densities  approximately  that  of  ion 
densities  exiting  thruster;  suggests  facility  pressure  remains  above  optimal 

*  Future  Work:  see  next  page 

—  Need  to  understand  the  dynamics  within  the  critical  regions  of  a  Hall  effect 
thruster  discharge  and  plume... 
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Critical  Regions  of  a  Hall  Effect 
Thruster  Discharge  and  Plume 


Division  of  Hall  effect  plasma  into  three  overlapping  regions 


Much  prior  work  has  examined 
one  of  more  aspects  of  facility 
effects  on  Hall  effect  thruster 
testing,  most  focused  on 
empirical  mitigation  of  effects 


Systematic  studies  of  the 
effects  of  ground  test 
facilities  on  each  of  these 
regions  is  necessary 


Discharge  (within  green  box) 

—  Region  of  ionization  and  acceleration  of  propellant 

—  Ranges  from  the  anode  to  vicinity  of  cathode 


Mid-Plume  (within  white  box) 

—  Luminous,  collisional  region  outside  thruster  body,  includes  near-plume  external  acceleration 

—  Region  with  isolated  slow  ion  population  in  central  core  and  responsible  for  cathode  coupling 


Far-Plume  (outside  yellow  trapezoid) 

—  Region  of  low  collisional  probability,  far  from  thruster  and  influence  of  applied  fields 

—  Region  of  interest  to  spacecraft  integrators  and  dependant  on  behavior  of  both  prior  regions 
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